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I.  Objectives 


This  is  the  final  perfonnance  report  for  our  grant  entitled  “Development  and  Simulation 
Studies  of  a  Novel  Electromagnetics  Code”.  The  main  purpose  of  this  research  program  is 
to  develop  novel  simulation  tools  for  the  design  and  study  of  high-power  microwave 
sources.  One  of  the  most  challenging  aspects  of  high-power  microwave  source  physics  is 
understanding  and  correctly  modeling  the  generation  and  subsequent  effects  of  space- 
charge  fields  in  the  source.  The  accurate  modeling  of  these  fields  is  one  of  the  key 
components  of  this  research  program.  Moreover,  the  subject  of  high-power  microwave 
sources,  as  a  whole,  has  been  of  great  interest  and  value  to  the  Air  Force  for  many  years. 
Hence,  this  research  program  is  of  great  relevance  for  the  Air  Force. 

The  research  program  has  four  facets: 

1.  Develop  novel  theoretical/numerical  tools  for  modeling  electromagnetic  fields 
due  to  space-charge  dominated  electron  beams  in  high-power  microwave  sources 

2.  Utilization  of  high-speed  computational  methods  to  maximize  the  computational 
efficiency  of  these  algorithms 

3.  Benchmarking  of  these  codes  to  determine  the  accuracy  of  the  numerical  methods 
as  well  as  to  identify  the  primary  errors  which  are  present  in  the  methods 

4.  Training  of  graduate  students  in  the  techniques  of  theoretical/computational 
electromagnetics 


II.  Status  of  Effort 


Within  the  past  three  years,  we  have  made  significant  progress  in  our 
development/implementation/benchmarking  of  novel  electromagnetic  algorithms  for  use  in 
high-power  radiation  source  simulations.  Specifically,  we  have  successfully  developed  a 
state-of-the-art  space-charge  solver  that  can  be  utilized  to  model  space-charge  fields  due 
to  high-current  electron  beams  in  either  cylindrical  or  rectangular  waveguides.  Such  a 
fonnalism  is  ideally  suited  for  modeling  the  drift  sections  of  klystrons  and  the  interaction 
regions  of  free  electron  lasers.  This  fonnalism  is  based  on  a  time-dependent  Green’s 
function  method,  and  has  the  primary  advantage  that  it  is  dispersion  free.  Dispersion 
errors  are  commonly  found  in  grid  based  electromagnetic  solvers,  such  as  the  Yee 
algorithm,  and  can  contribute  to  unwanted  errors  in  simulations,  such  as  numerical 
Cherenkov  radiation.  The  reason  why  our  method  is  dispersion  free  and  does  not  suffer 
from  these  unwanted  errors  is  because  the  Green’s  functions  are  constructed  from 
solutions  to  the  Helmholtz  equation.  Solutions  of  the  Helmholtz  wave  equation  have 
exactly  the  correct  dispersion  properties.  Hence,  any  linear  combination  of  these 
solutions,  such  as  the  Green’s  function,  also  has  the  correct  dispersion  properties.  An 
additional  advantage  of  our  method  is  that  since  the  Green’s  function  is  generated  from 
delta  function  sources  in  time  and  space  in  the  wave  equation,  it  can  be  readily  used  to 
model  the  space-charge  fields  of  arbitrarily  small  beam  bunches.  This  is  ideal  for  modeling 
small  electron  bunch  lengths  in  klystrons,  or  even  microbunching  instabilities  which  are 
commonly  found  in  free-electron  lasers. 

The  two  main  challenging  aspects  of  this  work  was  the  theoretical  development  of  the 
Green’s  function  method  and  the  computational  requirements  for  running  Green’s 
function  based  algorithms.  In  the  course  of  this  research  program,  we  have  discovered 
that  the  Green’s  function  method,  while  extremely  powerful  for  certain  simplified 
conductor  geometries,  has  the  limitation  that  it  is  not  apparent  how  to  extend  these 
methods  to  arbitrary  conductor  geometries.  During  this  research  program,  we  were  able 
to  develop  and  implement  Green’s  function  based  schemes  for  circular  and  rectangular 
conductor  pipe  geometries.  This  has  led  to  a  highly  accurate  description  of  space-charge 
effects  in  these  geometries,  which  has  applications  to  a  variety  of  high-power  microwave 
source  problems,  such  as  beam  dynamics  in  drift  sections  of  klystrons.  However,  the 
problem  of  how  to  extend  Green’s  function  methods  to  conductor  geometries,  such  as  the 
A-6  magnetron,  is  still  an  open  question.  The  computational  requirements  for  this 
scheme  have  also  shown  to  be  quite  challenging.  Unlike  in  a  typical  grid  based  scheme 
which  only  utilizes  the  values  of  electric  and  magnetic  fields,  as  well  as  charge  density  and 
current  density,  from  the  previous  time  step,  the  Green’s  function  method  uses  the  values 
of  fields  and  charge  and  current  densities  from  all  previous  points  in  time.  These 
requirements  lead  to  significant  challenges  on  detennining  the  optimum  algorithm  design. 

Our  significant  accomplishments  have  led  to  the  publication  of  a  paper  in  IEEE 
Transactions  on  Plasma  Science  entitled  “Exact  Space-Charge  Field  Solutions  for 
Cylindrically  Symmetric  Beam  Currents  in  a  Circular  Conductor  Pipe”,  the  training  of 
five  graduate  students  in  computational  electromagnetics,  and  the  completion  of  a  Ph.D. 


thesis  of  one  graduate  student,  Chong  Shik  Park.  Our  specific  scientific  achievements 
include  the  development  of  space-charge  field  solvers  for  cylindrical  and  rectangular 
conductor  geometries,  and  the  benchmarking  of  these  schemes.  We  have  also  developed  a 
parallelized  version  of  the  cylindrical  field  solver  method  using  MPI  (Message  Passing 
Interface)  scheme.  This  enables  our  Green’s  function  based  algorithms  to  operate  over  a 
large  number  of  simultaneous  processors.  Finally,  we  have  also  performed  a  novel 
calculation  and  initial  numerical  study  of  a  dispersion  free  FEL  (free-electron  laser)  model 
for  beams  in  a  rectangular  waveguide.  We  believe  that  this  result  may  be  of  great  interest 
to  the  Air  Force  for  future  FEL  applications. 


III.  Accomplishments/New  Findings/Future  Work 


Throughout  the  course  of  our  research  program,  we  have  made  a  number  of  important 
accomplishments  regarding  the  development  and  implementation  of  novel  space-charge 
field  solvers.  Our  first  major  accomplishment  was  the  theoretical  development  of  the 
Green’s  function  based  algorithm.  This  work  was  largely  outlined  in  our  IEEE 
Transactions  on  Plasma  Science  publication  in  2010,  as  well  as  Chong  Shik  Park’s  Ph.D. 
thesis.  For  a  cylindrical  conductor  geometry,  we  were  able  to  calculate  the  exact 
electromagnetic  space-charge  fields  given  an  arbitrary  cylindrically  symmetric  charge  and 
current  density,  which  satisfies  the  continuity  equation.  This  work  is  an  enormous 
theoretical  feat,  since  it  required  the  derivation  of  space-charge  fields  from  arbitrary 
sources,  which  is  a  highly  nonlinear  problem.  However,  the  beauty  of  this  scheme  is  that 
by  writing  down  the  solution  one  can  immediately  implement  it  numerically  to  calculate 
the  space-charge  fields  for  a  given  charge  and  current  density. 

We  utilized  this  scheme  to  investigate  the  space-charge  fields  of  a  bunched  beam  in  a 
pipe,  which  was  undergoing  radial  oscillations.  These  oscillations  could  be  caused  by  a 
TM0|  mode  in  the  1  GHz  range  having  a  characteristic  transverse  radial  electric  field  of  14 
MV/m,  which  is  interacting  with  the  beam.  This  calculation  served  two  purposes,  one  it 
allowed  us  to  explore  the  importance  of  fully  electromagnetic  space-charge  field  solvers 
versus  electrostatic  solvers,  as  well  as  to  act  as  a  benchmark  between  the  electrostatic 
method  and  our  electromagnetic  Green’s  function  method. 


Figure  1 


Figure  1  shows  plots  of  the  normalized  radial  electric  field  versus  a  normalized  radius 
for  the  Green’s  function  electromagnetic  method  (dotted)  and  a  Green’s  function  based 
electrostatic  method  (solid),  at  different  points  during  the  radial  beam  oscillation,  t=0, 
t=0.25T,  t=0.50T,  and  t=0.75T,  where  T  is  the  period  of  oscillation.  Details  of  the 
electrostatic  model  and  the  numerical  calculation  are  given  in  our  IEEE  paper,  as  well  as 
C.S.  Park’s  Ph.D.  thesis.  We  found  for  these  parameters  that  the  difference  in 
electrostatic  and  electromagnetic  space-charge  methods  can  vary  as  much  as  10%.  This 
may  be  significant  for  problems  involving  space-charge  dominated  beam  transport 
including  high-current  klystron  beams  in  a  drift  section. 

As  mentioned  before,  the  Green’s  function  method  has  a  significant  challenge 
compared  to  other  schemes,  such  as  the  Yee  algorithm,  in  that  it  requires  the  calculation  of 
integrals  over  charge  and  current  densities  for  all  previous  time  steps.  This  implies  that 
for  a  large  number  of  simulation  time  steps,  Nx,  the  total  computational  time  scales  as 
Nx.  In  order  to  improve  the  computational  efficiency  of  this  scheme,  we  developed  a 
Green’s  function  based  algorithm  which  utilizes  the  Message  Passing  Interface  (MPI) 
scheme. 

The  MPI  scheme  enables  one  to  distribute  the  computational  problem  over  many 
processors.  The  way  we  accomplished  this  is  by  implementing  a  grid  into  our  code.  The 
only  purpose  for  this  grid  is  to  keep  an  inventory  of  the  charge  and  current  densities  for 
all  time.  Since  we  are  not  solving  the  space-charge  fields  using  a  gridded  algorithm,  our 
code  still  maintains  its  dispersion  free  property.  Each  processor  is  “responsible”  for  a 
certain  range  of  cells  in  the  longitudinal  direction.  Specifically,  each  processor  performs 
integrations  over  the  charge  and  current  densities  in  its  range  of  cells.  Figure  2  shows  a 
plot  of  the  speed  up  of  the  algorithm  using  the  MPI  scheme  as  a  function  of  the  number 
of  processors.  One  can  immediately  see  that  we  can  achieve  excellent  performance  with 
the  MPI  scheme,  and  the  effects  of  crosstalk  between  processors  does  not  seem  to  be  a 
major  limitation  for  our  algorithm. 


Speed  Up  Curve(20,000  time  steps) 


Figure  2 


In  addition,  we  have  also  recently  tested  the  usage  of  Green’s  function  based  space- 
charge  field  solvers  for  modeling  free  electron  laser  systems  within  a  rectangular 
waveguide  structure.  Specifically,  we  have  developed  a  model  for  electron  beam  motion 
through  an  undulator  magnet,  and  used  our  dispersion  free  methodology  to  calculate  the 
space-charge  fields  within  the  beam.  Figure  3  shows  a  simulation  plot  of  the  space-charge 
electric  field  in  the  x-direction  for  a  highly  relativistic  electron  beam  moving  in  the  z- 
direction  with  an  energy  of  17.5  GeV,  a  bunch  charge  of  1.0  nC,  a  bunch  length  of  25pm, 
and  a  transverse  emittance  of  1.4  mm-mrad.  The  undulator  magnet  has  a  maximum 
magnetic  field  of  1.0  T  and  an  undulator  period  of  35.6  mm.  Although  these  initial  code 
test  parameters  are  associated  with  the  design  of  a  future  X-ray  FEL,  such  as  the 
European  XFEL  experiment,  it  is  straightforward  to  modify  the  parameters  in  our 
formalism  to  model  millimeter  wavelengths  that  are  more  pertinent  to  the  Air  Force.  The 
space-charge  field  is  plotted  in  the  beam  rest  frame  as  a  function  of  normalized  x  and  z. 
Utilizing  our  dispersion  free  scheme,  we  are  able  to  obtain  highly  accurate  space-charge 
field  calculations  for  FEL  systems.  This  may  prove  very  useful  for  future  Air  Force  free 
electron  laser  design  studies. 


Figure  3 


1.4 


Finally,  our  research  group  has  successfully  utilized  this  project  for  the  training  of 
graduate  students  at  Indiana  University.  This  project  has  supported  Chong  Shik  Park, 
who  recently  received  his  Ph.D.  in  2009,  Vivek  Krishna,  who  received  his  M.S.  degree  in 
computer  science,  Abhijeet  Kodgire  who  was  a  computer  science  student,  Adarsh  Yoga, 
who  received  his  M.S.  in  computer  science,  and  Rohith  Goparaju,  who  received  his  M.S. 


in  computer  science.  Dr.  Park  was  heavily  involved  in  the  development  of  the  dispersion 
free  theoretical  methods  and  numerical  implementation  of  these  methods.  Mr.  Krishna 
was  primarily  involved  in  developing  the  grid-based  algorithm  for  tracking  beam  charge 
and  current  densities  while  simultaneously  utilizing  our  dispersion  free  methodology  for 
computing  the  space-charge  fields.  A  short  time  after  Dr.  Park  and  Mr.  Krishna 
graduated  we  again  utilized  this  grant  to  begin  supporting  two  exceptional  computer 
science  graduate  students,  namely  Adarsh  Yoga  and  Abhijeet  Kodgire.  Mr.  Yoga  was 
responsible  for  benchmarking  our  dispersion  free  methods  against  a  finite-difference-time- 
domain  (Yee)  method.  Mr.  Kodgire  was  responsible  for  implementing  the  MPI 
framework  within  the  existing  grid-based  algorithm.  Mr.  Kodgire  left  our  group  in 
September  2010,  and  in  January  2011,  we  began  supporting  Rohith  Goparaju.  Mr. 
Goparaju  was  responsible  for  completing  the  MPI  implementation  of  our  code. 


IV.  Personnel  Supported 


Each  year  this  grant  has  included  funding  for  one  faculty  member  (M.Hess)  for  two 
months  summer  salary,  as  well  as  two  graduate  students,  with  stipend  and  fee  remission. 
One  graduate  student  supported  by  this  grant,  Chong  Shik  Park,  has  recently  graduated 
with  a  Ph.D.  in  physics  at  Indiana  University.  Dr.  Park  was  funded  at  the  start  of  phase  1 
on  this  grant  until  May  2010.  Another  graduate  student  supported  by  this  grant,  Vivek 
Krishna,  has  graduated  with  a  M.  S.  degree  in  computer  science.  Mr.  Krishna  had  been 
funded  since  the  start  date  of  Phase  2  of  this  grant,  i.e.  Jan.  1,  2009,  and  his  support  ended 
on  May  2010.  In  May  2010,  two  additional  graduate  students  in  the  computer  science 
department  were  supported  by  this  grant,  Adarsh  Yoga  and  Abhijeet  Kodgire.  In 
September  2010,  Mr.  Kodgire’s  support  ended  on  this  grant.  In  January  2011,  another 
computer  science  student,  Rohith  Goparaju,  was  added  to  this  grant.  Mr.  Yoga  and  Mr. 
Goparaju’s  support  from  this  grant  ended  on  May  2011. 


V.  Publications 


Peer  Reviewed  Journals 

M.  Hess  and  C.  S.  Park,  “Exact  Space-Charge  Field  Solutions  for  Cylindrically 
Symmetric  Beam  Currents  in  a  Circular  Conductor  Pipe”,  IEEE  Trans.  Plasma  Sci.,  38, 
pp.  1584-1591,2010 


Thesis 

C.S.  Park,  “Investigations  of  Electromagnetic  Space-Charge  Effects  in  Beam  Physics”, 
Indiana  University,  Ph.D.  Thesis  2009. 


VI.  Interactions/Transitions 


a. 


Talks: 

M.  Hess  and  C.S.  Park,  “Three-dimensional  Space-Charge  Modeling  of  Intense  Beams  in 
a  Circular  Conducting  Pipe”,  APS  Division  of  Plasma  Physics  Meeting  2009. 

M.  Hess  and  C.S.  Park,  “Simulating  Space-Charge  Physics  of  High-Current  Beams  with  a 
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Poster: 

M.  Hess,  A.  Yoga,  R.  Goparaju,  V.  Krishna,  A.  Kodgire,  and  C.S.  Park,  “Development 
and  Benchmarking  of  a  Novel  Grid  Based  Space-Charge  Field  Solver”,  International 
Conference  on  Plasma  Science  2011. 

M.  Hess,  C.S.  Park,  V.  Krishna,  and  A.  Kodgire,  "A  Dispersion  Free  Method  for 
Modeling  Space-Charge  Physics  in  a  Circular  Pipe",  International  Conference  on  Plasma 
Science  2010. 

M.  Hess,  “A  Dispersion  Free  Methodology  for  Modeling  Intense  Charged  Particle 
Beams”,  APS  Division  of  Plasma  Physics  Meeting  2010. 

M.  Hess  and  C.S.  Park  “A  Dispersion  Free  Three-Dimensional  Space-Charge  Modeling 
Method”,  Particle  Accelerator  Conference  2009. 


b.  There  were  no  consultative  or  advisory  functions  to  other  laboratories  or  agencies 
performed. 

c.  There  were  no  technology  assists,  transitions,  or  transfers  performed. 


VII.  New  Discoveries  /  Inventions  /  Patent  Disclosures 

None. 
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None. 
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Investigations  of  Electromagnetic 
Space-Charge  Effects  in  Beam  Physics 


In  order  to  achieve  high  quality  charged  particle  beams  for  electron  sources,  such 
as  rf  photoinjectors,  and  high-power  microwave  sources,  such  as  klystrons,  the  effects 
of  space-charge  need  to  be  addressed.  In  these  systems,  the  beam  energy  is  relatively 
low  while  the  beam  charge  is  relatively  high  leading  to  important  space-charge  physics. 
While  there  have  been  a  number  of  previous  works  which  have  analyzed  the  physics 
of  space-charge,  most  of  these  studies  have  looked  at  the  problem  assuming  that 
the  space-charge  physics  was  either  purely  electrostatic  so  that  analytical  techniques 
could  be  applied  or  fully  electromagnetic  using  only  computational  methods,  such 
as  the  Yee-PIC  algorithm,  for  simulating  these  effects.  The  present  study  is  unique 
in  that  it  incorporates  novel  theoretical  techniques  using  time-dependent  Green’s 
functions  for  computing  the  fields  and  modeling  the  space-charge  physics  analytically 
and  numerically. 

This  work  presents  the  calculations  of  time-dependent  electromagnetic  space- 
charge  fields  for  a  perfectly  conducting  pipe  with  and  without  a  cathode.  Using 
a  Lorenz  gauge  formalism  within  Maxwell’s  equations,  the  beam  space-charge  fields 
are  computed  from  a  time-dependent  Green’s  function  method.  In  this  method,  the 
correct  conductor  boundary  conditions  are  implemented  such  that  the  effects  of  image 
charges  and  image  currents  clue  to  the  cathode  and  cavity  walls  are  included. 

In  order  to  simulate  the  beam  dynamics  and  electromagnetic  space-charge  effects 
in  the  rf  photocathode  gun,  a  new  code  called  IRPSS  (Indiana  Rf  Photocathode 
Source  Simulator),  has  been  developed.  The  results  of  IRPSS  are  compared  with 
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electrostatic  codes,  such  as  a  time-independent  Green’s  function  based  space-charge 
algorithm,  as  well  as  known  beam  physics  codes,  such  as  PARMELA.  We  highlight 
the  important  numerical  challenges  and  computational  limits  which  we  have  analyzed 
during  the  development  of  IRPSS.  In  addition,  we  compare  IRPSS  results  to  an  ex¬ 
perimental  beam  loss  measurement  performed  on  the  Argonne  Wakefield  Accelerator 
rf  photocathode  gun. 

We  also  show  how  the  code  can  be  extended  to  include  important  physics  for 
beams  with  non- negligible  transverse  currents,  such  as  those  found  in  high-power 
microwave  sources. 
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Chapter  1 
Introduction 


In  most  of  electron  accelerator  systems,  such  as  FELs  (Free  Electron  Laser)  [1,  2] 
and  Colliders  [3,  4],  the  high  “quality”  electron  beam  has  been  critically  required. 
For  example,  the  FEL  needs  an  electron  beam  with  a  short  wavelength  and  dense 
phasespace  density  to  generate  coherent  radiation  sources.  In  the  high  energy  physics 
area,  the  e+e-  colliders  demand  a  high  number  of  possible  events  at  IPs  (Interaction 
Points).  In  these  applications,  it  is  very  important  to  achieve  a  bunched  beam  to  a 
small  spot  size,  i.e. ,  high  quality  beam. 

The  beam  quality  can  be  described  in  various  ways,  however,  the  beam  “bright¬ 
ness”  is  a  figure  of  merit  to  measure  the  beam  quality  for  the  electron  beam.  Qual¬ 
itatively,  the  brightness  is  a  measure  of  the  ability  of  a  beam  to  be  focused  in  6- 
dimensional  phasespace.  The  brightness  of  the  beam,  B,  is  defined  as  the  density  of 
particles  in  the  6-D  phasespace  [5]: 

^  spectral  flux  of  the  beam  2 Ipeak  ^  _ 

transverse  phasespace  area  en,xen,y' 

where  Ipeak  is  the  peak  current  of  the  beam,  and  en,x  and  en^y  is  the  normalized 
transverse  emittances.  Eqn.  (1.1)  tells  us  that  the  larger  peak  current  and  the  lower 
omittances  of  the  beam,  the  higher  the  beam  brightness  could  be  achieved. 
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The  emittance  in  Eqn.  (1.1)  is  the  area  occupied  by  the  beam  distribution  in  each 
transverse  phasespace  plane.  These  show  how  the  beam  is  focused  in  that  plane.  The 
term,  normalized,  includes  the  information  of  the  beam’s  energy.  The  normalized 
beam  emittance,  for  example  in  the  horizontal  direction,  is  given  by 

en,x  =  PlG  (1.2) 

where  /3  =  v/c,  7  is  the  relativistic  Lorentz  factor.  Eqn.  (1.3)  shows  the  normalized 
rms  emittance.  which  is  used  more  often  since  it  can  provide  the  information  of  the 
beam  regardless  of  the  particle’s  distribution  [6]. 

Cx  =  V(x2)ix'2)  ~  ixx')2-  (!-3) 

The  FEL  is  a  coherent  radiation  source  driven  by  an  electron  beam.  In  the  FEL, 
the  electron  beam  traverses  a  series  of  alternating  magnets,  such  as  a  wiggler  or  un- 
dulator,  and  radiates  light  at  wavelengths  depending  on  the  electron’s  energy,  wiggler 
period,  and  magnetic  field.  In  order  to  get  coherent  radiation,  the  electron  bunch 
must  be  focused  to  a  small  spot  size  less  than  or  equal  to  the  radiation  wavelength. 
This  requires  that  the  transverse  emittances  of  the  beam  must  be  small.  For  exam¬ 
ple,  SASE-FEL  (Self  Amplified  Spontaneous  Emission  FEL)  demands  the  normalized 
emittance  on  the  order  of  2  mm-mrad,  and  the  requirement  for  the  brightness  is 
9  x  1013  A/m2rad 2  [7].  The  4th  generation  of  X-ray  light  sources  needs  much  smaller 
emittances  (<  0.1  mm-mrad )  and  much  higher  brightness  (>  1017  A/m2rad 2)  [8]. 

In  high  energy  and  nuclear  physics,  accelerators  have  been  very  important  tools 
for  verifying  their  theories  and  models.  For  example,  the  Standard  Model  has  been 
successfully  proved  by  various  experiments  using  particle  colliders.  There  are  com¬ 
prehensive  efforts  to  find  the  Higgs  particle  to  complete  this  model.  Current  studies 
believe  that  it  requires  higher  than  100  GeV  of  center  of  mass  energy.  Recently,  there 
is  international  collaborative  efforts  to  design  and  construct  the  ILC  (International 
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Linear  Collider)  at  a  center  of  mass  energy  of  0.5  ~  1  TeV  [9].  This  machine  requires 
that  eUtX  ~  10  mm-mrad,  en)V  ~  0.02  mm-mrad,  and  B  ~  1013  A/m2rad2.  These 
parameters  are  very  challenging,  and  the  design  efforts  are  still  progressing. 

1.1  RF  Photoinjector 

Various  electron  sources  are  used  to  satisfy  the  requirements  for  parameters  of  ap¬ 
plications.  Usually,  the  rf  photoinjectors  are  widely  used  to  produce  high-brightness 
beams,  since  it  can  generate  a  bunched  beam  without  the  buncher  or  other  bunching 
devices  and  the  electron  beams  are  rapidly  accelerating  using  the  high  accelerating 
gradient  due  to  the  rf  field.  It  also  has  an  advantage  of  generating  an  initially  small 
omittance  beam  which  is  determined  by  the  wavelength  of  the  laser. 

The  rf  photoinjector  consists  of  three  main  parts,  the  photocathode,  the  laser 
system,  and  rf  cavities.  The  electron  beams  are  generated  by  hitting  the  laser  light 
to  the  photocathode,  and  the  they  are  accelerated  by  the  rf  electric  fields.  By  using 
the  high  gradient  rf  held,  the  beam  is  accelerated  more  rapidly  than  with  a  DC  held. 
The  hrst  rf  gun  was  tested  at  Los  Alamos  [10],  and  it  is  now  widely  used  in  numerous 
laboratories  as  electron  sources. 

The  rf  photoinjector  is  often  times  operated  in  a  TMoio  mode  of  a  pillbox  cavity, 
and  the  helds  in  each  cavity  are  out  of  phase  by  180  degrees,  i.e.,  7r-mode.  Figure  1.1 
shows  the  schematic  diagram  of  the  rf  photoinjector.  The  electron  beam  in  the  rf 
photoinjector  is  accelerated  under  the  influence  of  both  rf  and  space-charge  helds.  In 
general,  space-charge  helds  can  cause  significant  emittance  growth.  In  many  papers, 
people  have  tried  to  analyze  the  space-charge  effects  in  the  rf  photoinjector. 

K.-J.  Kim  hrst  studied  the  analytical  calculations  of  effects  due  to  rf  and  space- 
charge  helds  in  beam  dynamics  [11].  He  presented  that  (1)  the  acceleration  process 
due  to  the  rf  held,  (2)  the  effects  of  the  rf  acceleration  on  transverse  dynamics,  and 
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Figure  1.1:  Schematic  diagram  of  the  rf  photoinjector.  By  exciting 
the  photocathode  using  the  laser  (green),  the  electron  beam  (red) 
is  emitted.  The  beam  is  accelerated  by  the  rf  field  (black),  which 
is  operating  in  a  TMoio  7T-mode.  The  solenoid  magnetic  fields  are 
provided  to  reduce  the  emittance  growth  [12] 
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(3)  the  transverse  and  longitudinal  emittances  resulting  from  the  space-charge  effects. 
His  analyses  were  done  by  deriving  the  energy  scaling  of  the  field  components  under 
assumptions  of  a  quasi-electrostatic  approximation  in  the  rest  frame  of  the  beam  and 
a  I/72  dependence  of  the  space-charge  force. 

There  is  an  another  analytical  approach  to  investigate  the  space-charge  effects 
called  K-V  envelope  theory,  which  is  named  after  its  discovers,  Kaphchinskij  and 
Vladimirskij  [13].  Sacherer  calculated  the  space-charge  effects  in  the  accelerator  sys¬ 
tem  using  K-V  envelope  theory  [14],  This  theory  can  predict  the  size  of  the  beam  with 
both  space-charge  and  external  forces.  The  envelope  analysis  for  the  rf  photoinjector 
can  be  found  in  the  paper  of  Serafini  et  al  [15].  Although  it  can  provide  an  excellent 
description  of  the  beam,  however,  it  does  not  address  the  emittance  growth  resulting 
from  the  z-dependence  of  the  rf  and  space-charge  fields,  as  well  as  Wakefields,  due  to 
the  discontinuities  of  the  cavity  geometry. 

1.2  Main  Contributions  of  the  Thesis 

The  fully  analytical  treatments  of  the  electromagnetic  space-charge  fields  in  the  rf 
photoinjector  have  not  been  developed.  Since  the  beam  is  emitted  with  low  energy 
and  is  accelerated  to  a  relativistic  energy  in  a  few  cm,  its  beam  dynamics  are  strongly 
affected  by  the  space-charge  fields.  In  this  thesis,  our  goal  is  to  investigate  the  space- 
charge  effects  in  beam  systems.  Specifically,  we  will  focus  on  rf  photoinjectors  and 
high-power  microwave  sources.  The  analytical  and  numerical  solutions  of  electromag¬ 
netic  space-charge  fields  will  be  presented  in  the  presence  of  a  conducting  wall  and 
flat  cathode. 

Many  computational  codes,  such  as  PARMELA  [16],  ASTRA  [18],  TREDI  [17], 
and  MAFIA  [19],  are  capable  of  simulating  the  photoinjector.  However,  most  of  them 
utilize  Particlc-In-Cell  (PIC)  methods  to  calculate  the  fields  generated  by  the  beam. 
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In  order  to  understand  the  motion  of  the  beam  in  the  photoinjector,  we  need  to 
be  able  to  analyze  the  electromagnetic  fields,  since  the  space-charge  fields  strongly 
depend  on  time.  However,  the  PIC  algorithms  have  limitations  for  analyzing  the 
space-charge  fields.  And  there  are  unphysical  effects,  which  arise  in  PIC  codes,  for 
example,  the  numerical  grid  dispersion  and  numerical  Cherenkov  radiation  [20,  21]. 

Numerous  papers  have  studied  space-charge  effects  in  the  rf  photoinjector.  For 
example,  E.  Colby  et  al  showed  that  the  electrostatic  space-charge  induced  beam 
dynamics  for  a  10  psec  bunch  with  a  beam  radius  of  1  mm  disagrees  with  elec¬ 
tromagnetic  calculations  by  20%  at  locations  which  are  a  few  centimeters  from  the 
cathode  [22].  These  results  were  further  verified  by  A.  Candel  et  al  which  showed 
these  differences  and  found  that  the  discrepancy  can  be  larger  at  distances  beyond 
5  cm  for  the  LCLS  gun  [23].  Other  results,  such  as  by  J.-M.  Dolique  et  al,  suggest 
the  importance  of  the  self-generated  magnetic  space-charge  field  in  beam  dynamics 
within  the  rf  gun,  even  for  long  bunch  lengths,  i.e.,  100  psec  [24],  Near  the  cathode, 
one  expects  to  find  significant  differences  between  electrostatic  and  electromagnetic 
space-charge  solvers.  And  a  full  investigation  of  these  differences  near  the  cathode  is 
important  for  a  complete  physics  understanding  of  the  rf  gun. 

There  are  several  studies  which  analyzed  the  electromagnetic  fields  based  on 
Green’s  function  methods  [25,  26,  21].  These  references  [25,  26]  utilize  the  Fourier 
A;-space  in  the  longitudinal  direction  and  Bessel  function  expansion  in  the  transverse 
directions.  They  assumed  that  the  beam  has  a  uniform  density  and  is  accelerated 
uniformly.  But  the  uniform  density  beam  is  not  entirely  correct  during  the  acceler¬ 
ation  process,  and  the  beam  is  not  accelerated  uniformly  since  the  rf  held  is  varying 
in  time. 

In  an  effort  to  improve  the  modeling  capability  in  the  rf  photoinjector,  we  have  de¬ 
veloped  the  computational  code,  IRPSS  (Indiana  University  RF  Photoinjector  Source 
Simulator),  for  computing  the  electromagnetic  space-charge  fields  based  on  Green’s 
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function  methods.  This  code  calculates  the  electromagnetic  space-charge  potentials 
and  fields  numerically. 

In  application  for  high-power  microwave  sources,  an  extended  Green’s  function 
technique  will  be  discussed  in  Chapter  5.  This  method  will  include  the  effect  of 
transverse  currents  of  bunched  charged  particle  beams.  We  have  used  our  Green’s 
function  method  to  numerically  simulate  the  space-charge  field  of  a  bunched  electron 
beam  undergoing  “radial  breathing  mode  oscillations”  which  can  occur  in  a  high- 
power  klystron  drift  tube. 

1.3  Outline  of  the  Thesis 

Chapter  2  will  present  the  various  numerical  algorithms,  which  are  capable  of  model¬ 
ing  the  space-charge  fields  of  the  electron  beam  in  the  conducting  cavity,  especially, 
the  rf  photoinjector.  First,  we  will  show  the  existing  and  mostly  used  algorithms, 
such  as  SCHEFF,  Lienard-Wiechert  potentials,  and  Yee/PIC,  for  comparisons  with 
our  work.  As  mentioned  previously,  another  Green’s  function  approach  from  Salah 
et  al  will  be  discussed. 

Chapter  3  will  show  how  to  calculate  the  electromagnetic  space-charge  fields  in 
an  arbitrary  cross-section  conducting  pipe  with  longitudinal  beam  currents  in  the 
presence  of  the  cathode  at  z  =  0,  using  Green’s  function  methods.  For  specific  cases, 
we  will  present  examples  of  the  circular  pipe  cross-section  with  longitudinal  currents 
and  the  rectangular  pipe  cross-section  with  arbitrary  currents. 

For  the  case  of  the  circular  pipe  with  longitudinal  currents,  which  is  an  example  of 
the  rf  photoinjector,  the  simulation  code,  IRPSS,  has  been  developed.  In  Chapter  3, 
we  will  describe  the  computational  requirements  for  the  space-charge  field  calculations 
using  IRPSS.  The  benchmark  study  will  be  presented  to  check  the  validity  of  our 
Green’s  function  technique  and  numerical  implementations 
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In  Chapter  4,  we  apply  Green’s  function  techniques  to  compute  the  space-charge 
fields  in  the  rf  photoinjector.  The  parameters  of  the  ANL  AWA  1.3  GHz  rf  pho¬ 
toinjector  [27]  are  used  for  the  simulation  in  this  chapter.  We  will  discuss  how  the 
electromagnetic  space-charge  fields  near  the  cathode  are  different  from  the  electro¬ 
static  space-charge  fields  in  time.  With  these  space-charge  fields,  beam  dynamics 
simulations  are  conducted.  Specifically,  we  will  track  the  radial  phasespaces  for  both 
cases.  We  also  describe  the  results  of  an  experimental  beam  loss  measurement  which 
attempts  to  distinguish  the  fully  electromagnetic  space-charge  offsets  in  the  longitu¬ 
dinal  direction  from  the  electrostatic  space-charge  models. 

Chapter  5  presents  the  space-charge  field  solutions  in  the  presence  of  transverse 
currents  in  circular  conducting  pipe,  which  were  ignored  in  the  previous  chapter. 
The  electromagnetic  fields  are  derived  using  the  Bessel  function  series  expansions 
of  charge  and  current  densities  and  components  of  the  space-charge  fields.  These 
solutions  can  be  also  utilized  for  the  rf  photoinjector.  As  an  example,  we  will  see 
how  the  space-charge  fields  are  varying  for  a  radially  oscillating  beam  in  a  circular 
conducting  pipe. 

Finally,  in  Chapter  6,  we  will  conclude  this  dissertation  work.  And  we  will  discuss 
the  future  plan  for  improving  our  work. 
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Chapter  2 

Theoretical  Space-Charge 
Modeling 


There  have  been  numerous  codes  and  methods  developed  to  accurately  compute  the 
space-charge  effects  in  the  accelerator  system.  Especially,  in  order  to  design  the 
rf  photoinjector,  PARMELA  (Phase  And  Radial  Motion  in  Electron  Linear  Accel¬ 
erators)  [16]  is  a  workhorse  code  and  is  widely  used  for  various  examples.  It  can 
simulate  the  beam  dynamics  within  the  external  fields,  such  as  the  rf  and  solenoid 
fields.  However,  it  can  only  calculate  space-charge  fields  electrostatically.  TREDI 
(Three-Dimensional  Injectors)  [17]  is  fully  3-dimensional  simulation  code  for  charge 
particles  in  the  electron  source.  The  space-charge  fields  are  computed  using  Lienard- 
Wiechert  potentials.  PIC  (particle-in-cell)  is  the  one  of  the  most  commonly  used 
methods  to  simulate  the  beam  dynamics,  and  fully  electromagnetic  PIC  codes,  such 
as  MAFIA  [19],  compute  the  space-charge  fields  in  the  arbitrary  conducting  bound¬ 
aries  by  using  grid  based  field  solving  methods,  for  example,  Yee  algorithm  [28].  There 
are  also  analytical  approaches  to  compute  space-charge  fields  directly  from  Maxwell’s 
equations  using  either  Lorenz  or  Coulomb  gauge  condition  [25,  26].  In  this  chapter, 
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we  present  commonly  used  existing  space-charge  algorithms,  especially,  for  designing 
the  rf  photoinjector  system.  And  we  briefly  explain  our  space-charge  algorithm. 

2.1  Existing  Space-Charge  Algorithms 

2.1.1  SCHEFF  Algorithm  (PARMELA  and  ASTRA) 

PARMELA  and  ASTRA,  which  are  widely  used  codes  for  modeling  the  rf  photoin¬ 
jector,  use  an  electrostatic  based  space-charge  field  ring  algorithm  called  SCHEFF 
(Space  CHarge  EFFect)  [29].  In  PARMELA,  the  SCHEFF  algorithm  has  two  ver¬ 
sions,  SCHEFF  (2D),  and  SPCH3d.  In  this  present  work,  since  we  concentrate  on 
the  cylindrically  symmetric  system,  we  will  not  discuss  the  SPCH3d  code,  which  is 
a  3-D  PIC  routine,  and  calculates  the  field  solutions  based  on  a  convolution  of  the 
charge  density  and  the  Green’s  function  of  the  potential  [30].  There  are  also  ad¬ 
vanced  modifications  for  the  elliptically  symmetric  system  [31].  The  2-dimensional 
SCHEFF  algorithm  assumes  that  the  system  is  azimuthally  symmetric.  The  SCHEFF 
algorithm  breaks  up  the  macro-particles,  which  are  in  the  beam  volume,  into  a  set 
of  annular  rings,  and  integrates  trajectories  of  these  rings  through  the  fields  using 
the  phase  (time)  as  an  independent  variable.  Figure  2.1  shows  how  these  rings  are 
represented  in  the  grid  cells. 

In  order  to  include  relativistic  effects  of  accelerating  beams  in  the  SCHEFF  al¬ 
gorithm,  the  electrostatic  space-charge  fields  are  calculated  with  rings  in  the  average 
rest  system  of  the  bunch,  and  transformed  back  into  the  laboratory  frame,  again. 
The  electric  components  are  calculated  at  each  observation  point  by  performing  a 
two-dimensional  integration  in  (r,z  space),  where  each  element  of  integration  is  a 
circular  ring  of  charge  that  intersects  the  cell.  The  electric  held  components  at  the 
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Figure  2.1:  The  SCHEFF  method  models  the  electron  bunch  (a) 
as  a  collection  of  macro-particles,  and  (b)  breaks  into  concentric 
annular  rings.  Each  ring  is  assumed  to  contain  a  uniform  charge 
distribution,  (d)  The  ring  can  be  described  on  a  2-d  grid  where  each 
cell  is  the  intersection  of  a  ring  with  a  radial  plane  (constant  6)  [31] 
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rest  frame  can  be  expressed  as  [32], 


Er{r,z )  = 


47T€0  27t rb  r  [z2  +  (r  +  n)2]1/2 
r  r2  —  (r2  —  r2\  1 

r4)  - 


(2.1a) 


(2.1b) 


E  (r  z)  =  — _ - _ 4 rbzE[k) _ ^  (2.1b) 

47re0  27rr b  [~2  _j_  (r  _|_  Vb )2]1^2  [z2  +  (r  —  rb)2]  ’ 

where  rb  is  the  beam  radius,  Q  is  the  bunch  charge,  k  =  [4ryr / (^2  +  (r  +  r;,)2]1^2  and 
Jl  and  E  are  the  complete  elliptical  integrals  of  the  first  and  second  kind,  respectively. 
The  size  of  the  ring  is  scaled  with  the  variation  of  the  beam  size  to  match  the  amount 
of  the  charge  in  the  each  ring. 

However,  one  can  observe  limitations  in  PARMELA  as  the  SCHEFF  subroutine 
is  used.  First,  since  it  assumes  that  the  current  density  is  zero  in  the  beam  rest 
frame,  it  can  only  calculate  the  space-charge  fields  electrostatically,  i.e.,  not  includ¬ 
ing  time  retardation  effects.  In  PARMELA,  particles  are  pushed  by  the  external 
rf  and  solenoids  fields,  which  are  computed  from  SUPERFISH  and  POISSON  [33], 
respectively,  in  the  region  where  the  beam  travels.  It  takes  into  account  the  beam 
space-charge  fields  induced  by  the  real  charge,  as  well  as  the  image  charge  in  the  cath¬ 
ode.  A  second  limitation  is  that  even  though  it  can  also  include  interactions  between 
the  beam  and  the  cavity  walls,  i.e.,  the  wakefields,  these  fields  cannot  be  computed 
internally,  but  must  be  imported  externally.  A  third  limitation  is  that  the  radial  force 
generated  from  the  azimuthal  self-induced  magnetic  field  Bg  is  neglected.  Moreover, 
the  accuracy  in  PARMELA  are  in  question  for  low  energies  near  the  cathode  for  the 
high  brightness  regime  [31]. 
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2.1.2  Lienard-Wiechert  Potentials  (TREDI) 


TREDI  is  a  multi-purpose  macro-particle  3D  Monte  Carlo  code,  which  can  handle  the 
electron  beam  in  the  various  accelerator  systems,  such  as  rf  photoinjectors,  linacs,  and 
bending  magnets  [17].  The  beam  dynamics  in  TREDI  is  obtained  by  integrating  the 
Lorentz  force  law  by  means  of  a  4th  order  Runge-Kutta  algorithm.  The  external  fields 
from  rf  cavities,  solenoids,  or  bending  magnets  can  be  computed  analytically,  loaded 
from  external  hies,  or  evaluated  from  the  on-axis  fields  using  the  series  expansion. 

The  space-charge  fields  induced  by  the  beam  is  calculated  using  Lienard-Wiechert 
potentials  [34],  and  are  given  by 


2" 

1  h  x 

(h  —  (fj  x  /? 

i  (ft-'3) 

_  _L 

1  - 

(3 

_ L_\ _ / _ _|_ _ L _ . 

4«»  (l  — (5-ft)3|R|  4,r£o  (l-/3.ft)V|2 


(2.2a) 


and 


B  =  -n  x  E, 


(2.2b) 


where  R  =  |x-x'|,  n  =  R/  |R|,  and  these  fields  are  calculated  at  the  retarded  time, 


t'  =  t 


|R(f')l 


(2.3) 


c 

Similar  to  PARMELA  and  ASTRA,  TREDI  assumes  that  the  velocity  of  the 
source  particle  does  not  change  on  a  time  scale  comparable  to  the  retarded  time,  and 
the  contribution  of  acceleration  fields  is  negligible.  Moreover,  the  fields  are  calculated 
in  the  rest  frame,  and  then  are  transformed  to  the  laboratory  frame  using  the  Lorentz 
transformation: 


Elab  ^Ej-est, 


(2.4a) 


and 

Biab  =  P  x  Elab.  (2.4b) 

However,  these  fields  are  computed  in  the  presence  of  a  hat  cathode  and  in  free-space, 
i.e.,  there  is  no  other  conductor  boundaries,  such  as  cavity  walls. 
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2.1.3  Yee/PIC  Method  (MAFIA) 

MAFIA  utilizes  the  Yee/PIC  method  to  model  the  electromagnetic  space-charge  fields 
in  the  accelerator  system  [19].  The  Yee  algorithm  is  a  FDTD  (Finite-Difference 
Time-Domain)  method  to  solve  Maxwell’s  curl  equations,  and  it  has  been  widely 
used  to  calculate  the  electromagnetic  fields  for  numerous  problems,  such  as  antennas, 
microwave  circuits,  and  accelerator  systems  [28,  35]. 

In  order  to  solve  the  time-dependent  Maxwell’s  curl  equations, 


„  ^  <9B 

VxE=-ar 

(2.5a) 

and 

<9D 

V*H  =  J  +  w. 

(2.5b) 

for  an  isotropic  medium, 

first,  Eqn. (2.5)  are  written 

in  the  differential  forms  in 

Cartesian  coordinates  as: 

d  Bx  _  8Ey  d  Ez 
dt  dz  dy 

(2.6a) 

dBy  dEz  dEx 
dt  dx  dz 

(2.6b) 

dBz  dEx  dEy 
dt  dy  dx 

(2.6c) 

dDX  dHZ  dEy 

dt  x  dy  dz  ’ 

(2.6d) 

dDy  dHx  dHz 

dt  y  dz  dx  ’ 

(2.6e) 

and 

dDz  _  dHy  dHx 

dt  dx  dy 

(2.6f) 

In  the  Yee  algorithm,  the  spatial  coordinates  (■ i,j ,  k)  are  defined  as  (iAx,jAy,  kAz), 
where  Ax,  Ay,  and  Az  are  grid  separations  in  Cartesian  coordinates.  And  then,  each 
component  of  the  electromagnetic  fields  in  Eqn.  (2.6)  can  be  written  as 


Fn(i,j,  k )  =  F{iAx,jAy,kAz,nAt ), 


(2.7) 
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where  F  can  be  any  of  Ex,  Ey,  Ez,  Bx,  By,  and  Bz,  At  is  the  time  step,  and  n  is 
the  time  index.  The  differential  equations,  Eqn.  (2.6),  will  be  computed  using  central 
differential  approximations,  for  example, 

dFn(i,j,  k )  Fn{i  +  1/2,  j,  k )  -  Fn{i  -  1/2,  j,  k) 

dx  Ax 

and 

dFn(i,j ,  k)  _  Fn+1'2(i,j,  k)  —  Fw~1/2(i,  j,  k) 
dt  At 

Since  we  are  using  the  central  difference  approximations,  Eqn.  (2.8)  has  a  second- 

order  accuracy  in  space  and  time. 

The  field  components  of  E  and  B,  in  Eqn.  (2.6)  are  spatially  offset  from  one 
another,  as  seen  in  Figure  2.2.  In  Figure  2.2,  The  secondary  lattice  edges  (red) 
connect  centers  of  primary  lattices  (blue).  The  electric  field  components  are  constant 
along  the  edges  of  the  primary  lattice  and  are  parallel  to  the  primary  lattice,  while  the 
magnetic  field  components  are  constant  through  the  faces  of  the  primary  lattice.  The 
magnetic  field  components  are  constant  along  the  edges  of  the  secondary  lattice  and 
are  parallel  to  the  secondary  lattice,  while  the  electric  field  components  are  constant 
through  the  faces  of  the  secondary  lattice. 

In  addition,  E  and  B  are  also  offset  in  time  by  At/ 2  using  the  leap-frog  algo¬ 
rithm  [36].  These  offsets  leads  that  the  change  of  E  in  time  depends  on  the  change  of 
B  in  space,  and  vice  versa.  Therefore,  the  new  E  in  time  can  be  calculated  from  the 
old  E  in  time  and  the  spatial  derivatives  of  B  in  the  surrounding  lattice.  In  a  similar 
fashion,  the  new  B  can  be  obtained. 

In  addition  to  the  order  of  approximation  of  the  discretization  of  the  differential 
equations,  there  are  several  other  factors  that  impact  the  accuracy  of  the  finite  element 
modeling.  Since  dimensions  of  the  source  are  small,  grid  sizes  of  the  cell  must  be 
sufficiently  small  enough  to  resolve  the  wave  propagation  through  the  cell.  This 
requires  large  number  of  cells  and  consequently  demands  intensive  computational 


(2.8b) 
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Ex 

Figure  2.2:  Discretization  of  the  system  into  cubes  and  the  position 
of  field  components  on  the  grid. 


resources.  The  FDTD  method  also  suffers  from  numerical  dispersion.  The  numerical 
phase  velocity  of  the  electromagnetic  wave  induced  by  the  beam  can  be  slower  than 
the  physical  one.  This  implies  that  the  particles  can  be  traveling  faster  than  the 
electromagnetic  radiation  induced  by  themselves.  Consequently,  this  violates  the 
causality  condition  and  gives  rise  to  numerical  Cherenkov  radiation  [20]. 

2.2  Green’s  Function  Approaches 

2.2.1  Mode  Analysis  and  Series  Expansions 

There  have  been  analytical  approaches  to  solve  the  electromagnetic  fields  in  a  rf 
photoinjector  [25,  26],  which  are  based  on  the  Green’s  function  method.  These  ap¬ 
proaches  solve  the  Maxwell’s  equation  in  the  presence  of  the  boundary  conditions, 
such  as  a  cathode  and  the  cavity  walls,  using  two  different  gauge  conditions.  And  the 
cylindrical  symmetry  is  used  to  simplify  the  problem.  However,  the  reflection  waves 
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from  the  wall  are  excluded  since  the  system  assumes  the  cavity  radius  is  much  greater 
than  the  beam  radius. 

First,  the  scalar  and  vector  potentials  are  assumed  to  satisfy  the  “Lorenz  gauge” 
condition  [25], 

v'a4I=°-  <2-9> 

Then,  this  implies  the  wave  equations  for  potentials  as: 

(V  -  A  =  -/X0J,  (2.10a) 

and 

<2'10b> 

The  particular  solutions  of  Eqn.  (2.10a)  can  be  found  by  using  the  series  expansion 
of  A,  then  the  coefficients  of  the  expansion  are  calculated  from  the  Fourier  transfor¬ 
mation  in  normal  modes  and  the  Green’s  function  method.  And  the  general  solution 
is  easily  obtained  from  the  homogeneous  equation  in  Eqn.  2.10a.  One  can  get  the 
scalar  potential,  0,  using  the  Lorenz  gauge  condition,  in  Eqn.  (2.9). 

Another  approach  utilizes  the  “Coulomb  gauge”  condition  [26],  V  •  A  =  0.  This 
yields  the  wave  equations, 

V2</>  =  —  —  p,  (2.11a) 

eo 

and 

(v2-§)a  =  ^J  +  1|w)  (2.11b) 

Eqn.  (2.11)  is  solved  by  using  the  integral  transform  and  Green’s  function  method. 

These  two  techniques  have  advantages  that  it  can  be  used  to  calculate  the  space- 
charge  fields  to  arbitrary  accuracy  for  given  beam  charge  and  current  densities  when 
there  is  a  sufficient  amount  of  eigenmode  numbers  and  Fourier  modes. 
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2.2.2  Analytical  Method  with  Green’s  Function 

Similar  to  the  methods  developed  by  Salah  et  al  [25,  26],  we  have  developed  analytical 
solutions  to  Maxwell’s  equations  for  the  rf  photoinjector  in  the  presence  of  boundary 
conditions,  such  as  the  flat  cathode  and  conducting  cavity  walls.  The  arbitrary  longi¬ 
tudinal  current  densities  of  the  beam  are  assumed  to  be  in  axially  boundary  conditions 
with  an  arbitrary  cross-section.  We  find  the  electromagnetic  space-charge  potentials 
and  fields  using  the  time-dependent  Green’s  function  methods  in  Chapter  3. 

The  Green’s  function  methods  explicitly  include  the  causality  condition,  which  is 
important  in  the  beam  dynamics  of  the  bunch  head  when  the  beam  is  emitted  from 
the  cathode.  And  the  time-retardation  effect,  which  can  be  a  delayed  interaction 
between  the  head  and  tail  of  the  bunched  beam,  is  also  included.  Since  Green’s 
functions  are  solutions  to  the  wave  equations  with  the  delta  function  source,  it  is 
possible  to  accurately  compute  the  electromagnetic  fields  of  the  bunch  with  arbitrary 
bunch  length  in  the  longitudinal  direction.  In  addition,  since  the  space-charge  fields 
using  the  Green’s  function  methods  can  be  found  everywhere  in  space  and  time,  as 
opposed  to  a  grid  based  space-charge  algorithm,  this  method  has  an  advantage  of 
being  dispersion  free. 

In  Chapter  5,  we  will  show  the  more  general  solution  for  the  electromagnetic 
fields  in  the  pipe  structure  with  a  cylindrically  symmetric  beam  when  transverse 
currents  are  present.  Numerical  implementations  of  our  Green’s  function  approaches 
for  specific  examples  will  be  shown.  Finally,  in  order  to  check  the  importance  of  the 
electromagnetic  effects  in  both  cases,  we  also  derive  the  electrostatic  Green’s  function 
solutions.  These  will  be  discussed  briefly  in  the  following  chapters. 
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Chapter  3 

Space-Charge  Modeling  of  the  RF 
Photoinjector 

3.1  Electromagnetic  Space-Charge  Potentials  and 
Fields 

In  order  to  model  space-charge  physics  in  the  rf  photoinjector,  we  simplify  the  geom¬ 
etry  of  the  system  with  a  flat  cathode  which  is  located  at  z  =  0  and  a  surrounding 
conducting  pipe  of  an  arbitrary  uniform  cross  section  which  is  extended  to  infinity. 
An  electron  source  beam  which  is  generated  by  a  laser  pulse  consists  of  a  charge 
density  p(r,t)  and  a  current  density  J(r,  t).  Figure  3.1  shows  the  schematic  view  of 
this  system  with  the  electron  beam. 

In  this  section,  we  will  show  the  solutions  of  the  electromagnetic  space-charge 
potentials  and  fields  of  the  rf  photoinjector  system  with  an  electron  beam  for  two 
specific  cases.  First,  the  solutions  for  the  system  with  an  arbitrary  pipe  cross  section, 
but  with  only  longitudinal  beam  currents,  i.e.,  J(r,t)  =  Jz(r,t)ez,  will  be  shown.  As 
an  application  to  the  rf  photoinjector,  we  will  show  the  potentials  and  fields  of  the 
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Figure  3.1:  Schematic  of  an  electron  source  showing  a  planar  cath¬ 
ode  at  z  =  0,  a  surrounding  pipe  which  has  an  arbitrary  cross  section 
but  is  axially  uniform  for  z  >  0,  and  ac  accelerating  electron  beam 
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pipe  with  a  circularly  symmetric  cross  section.  Second,  we  will  present  the  solutions 
for  the  pipe  with  a  rectangular  cross  section  by  assuming  arbitrary  beam  currents. 
Although  the  system  with  a  rectangular  pipe  cross  section  is  only  occasionally  used 
in  electron  sources,  it  has  a  theoretical  interest  since  it  is  a  fully  solvable  system 
with  arbitrary  beam  currents.  The  electromagnetic  space-charge  fields  have  not  been 
solved  for  the  system  with  an  arbitrary,  uniform  pipe  cross  section  with  arbitrary 
beam  currents.  The  calculations  shown  in  this  chapter  have  been  outlined  in  our 
paper  [21] 

In  general,  the  beam  in  an  accelerating  structure  has  transverse  current  compo¬ 
nents,  as  well  as  a  longitudinal  current  component.  But  we  note  that  the  assumption 
of  the  beam  current  being  solely  in  the  longitudinal  direction  is  very  good  since  trans¬ 
verse  currents  will  be  typically  small  compared  to  the  longitudinal  current  in  the  rf 
photoinjector  system.  This  approximation  can  be  verified  by  comparing  the  ratio  of 
the  transverse  beam  velocity  to  the  longitudinal  beam  velocity.  When  the  beam  is 
emitted  from  the  cathode,  the  ratio  of  the  beam  velocity  can  be  approximated  by  the 
order  of  \v^_/vz\  ~  \E±/EZ\.  Each  electric  component  consists  of  both  the  external  rf 
held  and  space-charge  held  induced  by  the  beam  itself. 

For  the  ANL  1.3  GHz  photoinjector  [27],  the  longitudinal  accelerating  rf  held  is 
the  order  of  25  ~  100  MV jm.  The  transverse  rf  held  at  the  beam  edge  is  about 
3  %  of  the  longitudinal  electric  held.  Using  the  beam  parameters,  such  as  the  bunch 
charge  Q\,  —  1  nC,  the  beam  radius  r b  =  1  mm,  and  the  bunch  length  4  =  3  psec, 
the  transverse  space-charge  electric  held  is  less  than  10  %  of  the  rf  held.  As  the  beam 
is  accelerating  by  the  rf  held,  the  longitudinal  velocity,  vz,  is  reduced  by  a  factor  of 
l/y  due  to  the  relativistic  effect,  where  7  is  the  relativistic  Lorentz  factor.  Since 
the  velocity  ratio  is  less  than  10%  for  typical  photoinjector  systems,  it  is  sufficient 
to  approximate  the  beam  currents  for  electromagnetic  held  calculations  as  J(r,f)  = 
Jz(r,t)ez. 
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3.1.1  Arbitrary  Pipe  Cross-Section  with  Longitudinal  Cur¬ 
rents 


We  assume  that  the  system  has  the  cross  section  of  an  arbitrary  pipe  cross  section  with 
the  longitudinal  beam  currents,  i.e.,  the  current  density  is  given  by  J(r,  t)  =  Jz( r,  t)ez. 
For  a  given  charge  density  p(r,t)  and  longitudinal  current  density  J2(r,f),  the  elec¬ 
tromagnetic  potentials,  (f>(r,t)  and  A(r ,t)  =  Az(r,t)ez,  yield  the  electromagnetic 
space-charge  fields,  E(r,  t)  =  —  V0  —  dA/dt  and  B(r,  t)  =  V  x  A,  and  these  can  be 
solved,  i.e., 

( n2  p(r ,t) 


0(r’t)  = 


1  d2\ 

y2  -  ^2^2  )  Mr,t)  = 


with  the  Lorenz  gauge  condition, 

lf  +  V.A  =  °.  (3.3) 

The  beam  charge  and  current  densities  satisfy  the  continuity  equation, 

^  +  V  •  J  =  0.  (3.4) 

In  this  system,  the  flat  cathode  and  the  pipe  are  perfect  conducting  boundaries, 
and  therefore  the  fields  must  satisfy  the  boundary  conditions: 

E(E  *)  X  “(surface  =  B(r;  *)  ‘  “(surface  =  °>  (3‘5) 

where  h  denotes  the  normal  vector  on  the  conductor  surface.  Eqn.  (3.3)  can  be 
satisfied  by  specifying  the  boundary  conditions  for  the  potentials  to  be: 


A  y  ft\  —A  \z>  —  D 

^  A  ' 4 1  surface  —  I  surface  —  u’ 
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and 

f)  a  z=0 

(n  •  V)  (A  ■  n)|surface  =  — —  =0.  (3.8) 

^  surface 

We  note  that  Eqns.  (3.6),  and  (3.7)  yield  the  correct  boundary  conditions  for  the 
fields,  and  Eqn.  (3.8)  is  obtained  from  Eqns.  (3.6),  and  (3.7)  along  with  the  gauge 
condition  in  Eqn.  (3.3). 

The  electromagnetic  potentials  may  be  written  in  terms  of  two  electromagnetic 
Greens  functions  G $  (r,  t ;  r',  t ')  and  Gz  (r,  t;  r',  t')  as 

0  (r ,t)  =  —  [  f  Gfj,  (r,  t]  r',  t')  p  (r;,  t')  d3r'dt',  (3.9a) 

^0  J — oo  J 


0*^  =  0,  (3.11a) 

G..£Le  =  0.  (3-llb) 

and 

dGz  "=0 
c)z 

surface 

from  the  boundary  conditions  of  0  and  Az  in  Eqns.  (3.6),  (3.7),  and  (3.8).  Here, 
the  unprimed  coordinates  in  space  and  time  indicate  the  field  point,  while  the  primed 
coordinates  are  the  source  point. 

Since  the  Green’s  functions  which  are  defined  by  Eqn.  (3.10)  are  generated  by 
delta  functions  in  both  space  and  time,  it  is  possible  to  calculate  the  electromagnetic 


(3.11c) 
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potentials,  and  hence,  the  electromagnetic  fields  for  given  beam  charge  and  current 
densities  as  shown  in  Eqns.  (3.9a)  and  (3.9b). 

One  can  find  the  solutions  to  G^z  from  the  homogeneous  solutions  of  Eqn.  (3.10)  [37], 

i.e., 

r<P)  (f  —  +') 

G<p  =  - — Y  ^ run  (r±)  Vmn  (r±)  x  [r_  -  T+]  (3.12a) 

mn 

and 

rf)  (+  —  +')  . 

Gz  =  — — - Y^mn  (r±)  (r±)  x  [T-  +  r+] ,  (3.12b) 

mn 

where  rj_  =  r  —  zez  and  r'j_  =  r'  —  z'ez  are  transverse  vectors,  the  normalized 
eigenfunctions  (rj_)  are  the  Dirichlet  solutions  of  the  transverse  2D  Helmholtz 
equation  with  eigenvalues  /cj_mn,  i.e., 


(V  +  ^Imn)  (r_l)  —  0, 


(3.13a) 


and 


'tpmn  (rJ-)lsurf  —  0) 

| i>mn  (r±)\2  d2r±  =  1, 


gik(z±z  ) 


r+  = 


dk- 


CVklmn  +  k2  (t  —  if) 


71 


Vk. 


2  _i_  U2 

±mn  '  ™ 


(3.13b) 

(3.13c) 

(3.14) 


The  Green’s  functions  defined  by  the  eigenfunction  expansions  in  Eqns.  (3.12a) 
and  (3.12b),  as  well  as  the  associated  T_  and  T+  defined  by  Eqn.  (3.14),  can  be 
readily  obtained  in  the  following  manner.  First,  by  using  the  fact  that 


©  (t  -  t')  ipmn  (r_j_)  elkz  {Cmn  cos  [a;  (t  -  t')}  +  Smn  sin  [a;  (t  -  t')]} 


(3.15) 


forms  a  complete  set  of  eigenfunctions  which  are  solutions  to  the  homogeneous  3D 
Helmholtz  equation,  where  Cmn  and  Smn  are  constants  and  lu  =  c^/k\mn  +  k2,  we 
can  express  G as  a  sum  over  these  eigenfunctions.  The  factor  0  (t  —  t')  is  included 
to  ensure  causality.  By  inserting  the  expansion  for  G into  Eqn.  (3.10)  and  then 
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integrating  both  sides  by  dt  in  the  limit  e  — >  0,  it  can  be  easily  shown  that 
Cmn  =  0  and  G^z  and  Tj_  are  given  by  Eqns.  (3.12a),  (3.12b)  and  (3.14).  To  be 
specific,  the  terms  in  corresponding  to  T+  yield  an  extra  term  proportional  to 
8  (z  +  z')  on  the  right-hand  side  of  Eqn.  (3.10).  However,  since  z,  z'  >  0,  this  term  is 
zero.  The  factor  T+  in  Eqns.  (3.12a),  and  (3.12b)  is  necessary  to  satisfy  the  cathode 
boundary  conditions  at  z  =  0  given  by  Eqns.  (3.11a)  and  (3.11b).  From  a  physical 
point  of  view,  T_  and  T+,  represent  the  contributions  of  the  real  charge  and  image 
charge  (due  to  the  cathode  at  z  —  0),  respectively,  on  the  electromagnetic  fields. 
We  note  that  the  method  of  constructing  Green’s  functions  from  eigenfunctions  is  a 
well-known  method  that  is  used  extensively  for  solving  non-homogeneous  scalar  wave 
equations  [37]. 

The  integral  in  Eqn.  (3.14)  may  be  simplified  [38]  to 

r±  =  Jo  (k±mnX±)  0  (A|)  ,  (3.16) 

where  Jq(x)  is  the  zeroth  order  Besssel  function  of  the  first  kind,  Q(x)  is  the  Heaviside- 
step  function,  and 

A|  =  c2  (t  -  t'f  ~(z±  z'f  .  (3.17) 

Using  Eqn.  (3.4)  and  integration  by  parts,  one  can  readily  check  that  the  solutions 
given  by  Eqns  (3.9),  and  (3.12)  satisfy  the  Lorenz  gauge  condition.  We  note  that 
the  factors  0  (A2 )  and  0  (A+)  enforce  a  causality  condition  on  the  electromagnetic 
waves  emanating  from  the  beam  charge  and  induced  image  charges,  respectively. 
When  analyzing  the  potentials  near  the  front  of  the  bunch,  these  factors  allow  for 
rapid  numerical  convergence  since  only  sources  “sufficiently  close”  to  the  point  of 
observation  need  to  be  considered. 

The  electromagnetic  fields  E  (r,  t)  and  B  (r,  t)  can  be  computed  from  the  poten¬ 
tials  in  Eqns.  (3.9).  In  general,  the  electric  field  will  have  components  in  both  the 
longitudinal  and  transverse  directions,  but  under  the  assumption  of  only  longitudinal 
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beam  currents,  the  magnetic  field  will  only  have  components  in  the  transverse  direc¬ 
tions.  The  transverse  electric  and  magnetic  fields  follow  immediately  from  Eqns.  (3.9), 
and  are  given  by 

E±  =  Y  (r^)  f  J  ^ mn  (r±)  x  [T_  -  T+]  p  (r 't')  dh'dt',  (3.18) 


and 
Bj_  = 


~^Y^zX  W^rn.n  (r±)  J  J  ip'mn  (r'±)  x  [T_  +  T +]  J z  (r 't')  dh'dt'.  (3.19) 


The  final  solution  to  the  electric  field  in  the  longitudinal  direction  requires  a  few 
mathematical  steps.  To  start, 


Ez  =  — —  f  [  dh'dt' 
^0  J — oo  J 


dG,h  1  d  Gz 

P  H - -  »  d  7 


dz 


c2  dt 


x  [  dz' 


^  (r±)  (r±)  /  dt' 


dX _  <9T_  d\+  <9T+  \  (  dX _  <9T_  dX+  <9T+  \  Jz 

P+  I  + 


dt  dX-  dt  dX+ 


dt  <9A_  dt  dX+  I  c 2 


Then  using  the  fact  that  dJ0(x)/dx  =  —J\(x), 


dX± 

and  the  completeness  relation 


\z±z' 


(3.20) 


dB (Ab  =  2A±(S (A|)  =  A±<iy±c  +  c(i-o)]  +  y±(,-c(i-f-))]_  (3 21) 


Y  (r±)  VCn  (r±)  =  &  (r±  -  r±)  , 


mn 


(3.22) 
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we  can  express  the  longitudinal  electric  field  as 


Ez  ^  rJ_  ^  v  (r±)  ^mn  (r_l_)  j 


X  /  d ^  l  ,  ""  J\  (fclmnA-) 


A_ 


/  A  ,  ^  (t  -  *') 
-{z-z)p+- 


c 


0  (A2+) 


>1 1  (^’±mnA- 


-  (z  +  z')  p  + 


Jz  ( t  -  t') 


+~2 T  dt' 

Ze0  J- oo 


x  |0  [z  -  c  (t  -  t')}  [p  (r_i_,  z  —  c(t  —  t'),  t')  -  Jz  (rj_,  z  -  c  (t  -  t') ,  t')  / c] 

-  0  [z  +  c  (t  -  t')]  [p  (rj_,  2:  +  c  (t  —  t') ,  t')  +  Jz  (r_i_,  z  +  c(t  —  t') ,  t')  /c] 

—  0  [-2;  +  c  (t  -  t')}  [p  (r_L,  -z  +  c(t-  t ')  ,t')  +  Jz  (r_L,  —z  +  c(t  —  t')  ,t')  /c] 

(3.23) 


We  note  that  the  step  functions  in  the  last  integral  term  arise  from  the  fact  that 
the  delta  functions  in  Eqn.  (3.21)  will  only  yield  finite  contributions  after  integration 
in  z'  when  the  zeros  of  the  delta  functions  occur  for  z'  >  0.  Since  it  is  always  true 
that  —z  —  c(t  —  t')  <  0,  the  term  which  is  proportional  to  0  [—z  —  c  (t  —  t')}  does  not 
contribute  to  the  last  term  in  Eqn.  (3.23). 


3.1.2  Circular  Pipe  Cross-Section  with  Longitudinal  Cur¬ 
rents 

For  the  special  case  of  a  circular  pipe  with  radius  a,  the  transverse  eigenfunctions  and 
corresponding  eigenvalues,  which  solve  Eqns.  (3.13),  are  given  in  cylindrical  coordi¬ 
nates  by 

1  J m  (jmnr/a)  eime 

Oj\[t£  (j  ran)  | 


Ipmn  (r±) 


(3.24a) 
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and 


(3.24b) 


where  Jm(x)  is  the  mth  order  Bessel  function  and  jmn  is  the  nth  positive  root 
of  Jm{x).  Hence,  Eqns.  (3.24),  along  with  Eqns.  (3.18),  (3.19),  and  (3.23)  yield  a 
complete  description  for  computing  the  space-charge  fields  in  a  circular  pipe  with  a 
cathode  at  z  —  0. 

We  should  emphasize  that  cylindrical  symmetry  is  only  placed  on  the  pipe,  but 
not  on  the  charge  and  current  distributions.  Therefore,  it  is  possible  to  use  this 
formulation  to  simulate  the  electromagnetic  fields  for  arbitrary  longitudinal  beam 
currents,  such  as  those  with  dipole  and  quadrupole  moments,  within  a  cylindrical 
pipe. 

As  a  numerical  example,  we  choose  the  parabolic  beam  in  the  transverse  direction 
with  the  zero- thickness  in  the  longitudinal  direction.  The  exact  electromagnetic  fields 
are  given  by 


neo rla2  ^  jon\Jl(jon)\2  [jonh 


~J\  ( jonfb )  Jo  (jOn'f’fe) 


(3.25a) 


x 


dr' 


Jo  (ionA-  )  @(A_)  —  Jo  ( JOnA+  )  0(A+) 


U!  ^cQdO\^  J\  (jonf  ) 
Be  (r ,  z,  t)  =  >  .  ,J(. 

nrba  ^  JOn|Jl(jOn)h 


JO  ■rj'b 


'J\  (, jon.T’b )  Jo  (jOn'f’fe) 


X 


dr' 


Jo  ( .7  On  A—  )  0(A_)  +  Jq  ( JOnA+  )  0(A+) 


dz"  (r) 
dr 


(3.25b) 
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and 


Ex 


4 Q  sr-'  _Jo  ijonJ' ) 

n  2 


ne0rla2  ^  joJMjonW 


JO  nTb 


~J\  (,jonTb)  Jo  ( jonTb ) 


X 


X 


dr'  l  <  [z  —  z"(t')\  —  (t  —  r') 


dz"  (r) 
dr 


(jonA-)  0(A2_)  +  2a2J0  (jonA— )  <5(A2_) 
A_ 


—  <1  [z  +  z" {t')\  +  (t  -  r') 


x 


J  Or 

X 


d?"  (r) 
dr 


J\  ( Jon A-(- )  0(A+)  +  2a"  Jq  ( JonA+  )  <5(A+) 


(3.25c) 


In  Eqn.  (3.25),  the  spatial  and  temporal  coordinates  are  normalized  to  express  the 
data  of  the  gun  in  general  make  parameters  dimensionless.  The  coordinates  are 
normalized  by  the  radius  of  the  rf  cavity,  a,  such  that  f  —  r/a,  z  —  z/a,  and  t  =  ct/ a, 
where  c  is  the  speed  of  the  light.  And  f &  =  r^/a  is  the  beam  radius,  Q  is  the  bunch 
charge,  and  z"(r)  is  the  bunch  location.  The  detailed  derivation  of  these  fields  are  in 
the  Appendix  A. 

Figures  3.2  and  3.3  show  examples  of  held  plots.  The  space-charge  fields  are 
numerically  computed  with  the  parameters  of  BNL  2.856  GHz  gun.  Figure  3.2  shows 
a  3D  plot  of  Er  vs.  r  and  z.  There  are  two  distinguishing  features  in  Figure  3.2:  a 
large  pulse  near  z/a  =  0.225,  which  is  the  bunch’s  axial  location  and  a  second  smaller 
pulse  with  an  outer  edge  at  z/a  =  0.25.  The  pulse  at  z/a  =  0.225  represents  the  local 
electric  held  generated  by  the  beam,  while  the  second  pulse  represents  the  initial 
electromagnetic  waves  which  were  generated  at  t  —  0  when  the  beam  is  emitted. 
Since  the  time  in  Figure  3.2  is  t  —  0.25 a/c,  the  point  z/a  =  0.25  corresponds  to  the 
location  of  the  causality  condition,  z  =  ct.  From  Figure  3.2,  it  is  also  seen  that  the 
holds  are  zero  when  z  >  ct,  which  illustrates  the  effect  of  the  causality  condition. 

Figure  3.3  shows  a  plot  of  Ez  vs.  r  and  z.  Figure  3.3  shows  similar  features  as 
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E,(4*roa7/0) 


Figure  3.2:  Plots  of  Er  vs.  r  and  z,  for  t  =  0.25a/c. 
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Figure  3.3:  Plots  of  E~  vs.  r  and  z,  for  t  =  0.25 a/c. 
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Figure  3.2,  namely,  two  pulses  -  one  near  the  beam  and  one  generated  by  the  initial 
waves  due  to  bunch  emission,  i.e.,  the  transient  electromagnetic  pulse  clue  to  the 
transient  radiations.  The  sharp  discontinuity  in  the  pulse  near  the  beam  is  located 
precisely  at  the  bunch’s  location,  and  denotes  the  point  at  which  the  sign  of  the 
electric  field  changes. 


3.1.3  Rectangular  Pipe  Cross-Section  with  Arbitrary  Cur¬ 
rents 


In  contrast  to  the  previous  example,  the  explicit  solutions  of  the  electromagnetic 
space-charge  potentials  and  fields  to  the  rectangular  pipe  cross  section  with  the  ar¬ 
bitrary  beam  currents  can  be  easily  found.  In  this  system,  the  electron  beam  has 
transverse  currents  as  well  as  the  longitudinal  currents,  i.e.,  J(r,  t)  =  Jx(r,t)ex  + 
Jy(r,t)ey  +  Jz(r,t)ez. 

As  in  the  previous  case,  the  electromagnetic  potentials  0  (r ,t)  and  A(r,f),  may 
be  solved  in  the  Lorenz  gauge.  The  difference  in  this  case,  however,  is  that  A  (r,  t ) 
will  also  have  transverse  components  instead  of  only  the  longitudinal  component.  For 
p(r,f)  and  J(r,f)  which  satisfy  the  continuity  equation,  Eqn.  (3.4),  the  potentials 
are  given  by 


pM) 

eo 


(3.26a) 


and 


_1_ 

c 2  d t2  J 


A  fat) 


— //,,.I(r. /), 


(3.26b) 


and  satisfy  the  boundary  conditions  listed  in  Eqns.  (3.6),  (3.7),  and  (3.8). 

Similar  to  the  previous  case,  the  solutions  to  the  electromagnetic  potentials  can 
be  expressed  in  terms  of  four  time-dependent  Green’s  functions,  which  we  denote  by 
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Gi  (r ,t;  r where  i  —  <f),  x,  y,  and  z.  The  potentials  are  given  by 


(j)  (r,  t)  =  —  [  j  G $  (r,  t-  r',  t')  p  (r',  t')  d3r'dt\ 

eo  J- oo  J 


(3.27a) 


Aa  (r,  t)=  no  I  I  Ga  (r,  t;  r',  i')  Ja  (r',  t')  dVcft', 


(3.27b) 


where  a  —  x,  y,  and  z.  All  four  Green’s  functions  satisfy  the  differential  equation 


V2--—  Gj:  =  -5(r-r')5(t-t'), 


(3.28) 


From  the  boundary  conditions  as  in  Eqns.  (3.6),  (3.7),  and  (3.8),  the  Green’s 
functions  are  subject  to  the  following  conditions,  namely 


G<t>  I  surf 


(3.29a) 


TX\y=0,y=W  Qx 


(3.29b) 


I  x=0,  x=L 


I  _  ^ y  _  n 

^ V\x=0,x=L  a 

UiJ  y=0,y=W 


(3.29c) 


\x=0,x=L  _  q 

Uz \y=o,y=W  ~  ~  U’ 

z= 0 

as  we  specify  the  sides  of  the  pipe  by  x  =  0,  x  =  L,  y  —  0,  and  y  =  W. 


(3.29d) 


The  solutions  to  Gi  can  be  readily  constructed  in  a  fashion  similar  to  the  previous 


2c  .  /mi rx\  .  ( mnx'\  .  / niry\  .  fnny'\ 

G*  =  ~ew  Z  Z sin  (— ) sm  (— J sin  (if)  8111  yw ) |r- " r+l  ■ 

7n=l  n=  1  x  x  ' 


(3.30a) 


x  T.w  Z^  Z^ 


m= 1 n= 1 


/ rmrx\  l  mnx  \  .  (mry\  .  nny  \  r  n 

cos  t  cos  — r~  8111  ~nrr ) sm  ( ~nrr  )  -  r+]  > 


(3.30b) 


Gy  t,w  S1 


m=  1  7i—l 


.  / rrmx\  .  /  mnx  \  / nny\  l  nny  \  r_^  _  , 

sm  t  sin  “ T"  cos  ATT  cos  -FiH  -  r+l  > 


(3.30c) 
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Gz 


where 


k 


_L  mn 


2 


[r-  +  r+] . 

(3.30d) 

(3.31) 


The  reason  why  it  is  relatively  straightforward  to  construct  the  potential  solutions 
in  this  geometry  for  arbitrary  beam  currents  as  opposed  to  a  general  curvilinear  pipe, 
such  as  a  circular  pipe,  is  due  to  two  reasons.  First,  in  the  rectangular  case  for  which 
it  is  appropriate  to  utilize  Cartesian  coordinates,  Eqn.  (3.26b)  can  be  separated  into 
three  component  equations.  In  these  equations,  there  is  no  coupling  between  the  two 
transverse  directions,  i.e.,  Jx  generates  Ax  but  not  Ay,  and  vice  versa.  This  is  not 
the  case,  however,  in  general  curvilinear  coordinates,  in  which  both  Ar  and  Ag  are 
generated  when  either  Jr  or  Jg  are  present. 

Second,  in  order  to  satisfy  the  Lorenz  gauge  condition  with  the  potentials,  it  is  nec¬ 
essary  to  change  the  order  of  differentiation  in  the  eigenfunctions  from  unprimed  (field) 
coordinates  to  primed  (source)  coordinates  using  integration  by  parts,  and  then  to  ex¬ 
ploit  the  continuity  equation.  This  was  readily  accomplished  for  longitudinal  currents 
when  the  eigenfunctions  in  the  longitudinal  direction  and  in  time  are  sinusoidal  as 
in  Eqn.  (3.14).  Since  the  transverse  eigenfunctions  for  the  rectangular  case  are  also 
given  by  sinusoidal  functions,  the  gauge  condition  can  be  readily  satisfied  for  the 
transverse  potentials  as  well. 
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3.2  Numerical  Implementation 

3.2.1  Development  of  IRPSS  Code 

We  develop  the  computational  code,  IRPSS  (Indiana  University  Rf  Photoinjector 
Source  Simulator),  for  computing  the  electromagnetic  space-charge  fields.  IRPSS  is 
capable  of  solving  the  space-charge  fields  with  Green’s  function  methods  numerically. 

We  assume  that  the  system  is  cylindrically  symmetric,  so  that  we  can  use  the 
results  of  the  previous  section.  We  also  assume  that  the  beam  is  moving  through  the 
prescribed  trajectory,  which  is  calculated  with  the  presence  of  the  external  rf  field, 

ErJ  =  E'q  cos(kz)  sin  (cut  +  99),  (3.32) 

where  E0  is  the  peak  rf  electric  field,  u  is  the  rf  frequency,  and  ip  is  the  injection  rf 
phase  [39].  This  implies  that  there  is  only  the  longitudinal  beam  currents.  Hence¬ 
forth,  the  electromagnetic  space-charge  fields  are  computed  for  the  bunch  with  the 
prescribed  trajectory. 

The  numerical  examples  in  this  work  are  using  the  parameters  of  the  two  rf  pho¬ 
toinjectors,  which  are  the  BNL  2.856  GHz  and  the  ANL  AWA  1.3  GHz  guns.  The 
main  gun  parameters  are  given  in  the  Table  3.1. 

The  distribution  of  the  bunched  beam  can  be  generalized  to  several  types  of  dis¬ 
tribution  functions,  such  as  parabolic,  Gaussian,  and  waterbag  in  3-dimensional  di¬ 
rections.  In  particular,  the  charged  particle  beam  distribution  in  the  rf  photoinjector 
has  a  short  bunch  length,  which  is  much  smaller  than  the  transverse  coordinates  of 
the  beam.  First,  we  assume  that  the  beam  has  the  finite  transverse  size  with  the  zero 
thickness  in  the  longitudinal  direction,  and  then  the  bunched  beam  consists  of  many 
sliced  (zero  thickness)  beams.  The  multi-sliced  bunched  beam  will  be  discussed  in 
the  following  Section. 
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BNL 

ANL 

rf  frequency  (GHz) 

2.856 

1.3 

peak  rf  field  (MV /m) 

100 

25  ~  77 

rf  injection  phase  (deg) 

68 

65 

radius  of  the  cavity  (m) 

0.04111 

0.0908 

radius  of  the  beam  (mm) 

1 

1 

charge  of  the  beam  (nC) 

1 

1 

Table  3.1:  The  rf  photoinjector  beam  parameters  of  the  BNL  2.856 
GHz  and  the  ANL  AWA  1.3  GHz  guns 


The  space-charge  fields  are  determined  by  the  beam  profile,  which  must  be  up¬ 
dated  at  each  time  step  of  simulations.  In  this  work,  we  restrict  that  the  beam  density 
function  is  not  changed  with  time  and  uncoupled  between  the  longitudinal  and  trans¬ 
verse  coordinates,  then  it  can  be  separated  into  longitudinal-  and  time-dependent 
parts  and  a  transverse-dependent  one.  To  simplify  the  problem,  we  assume  that  there 
is  no  transverse  motion.  Using  the  completeness  of  the  transverse  eigenfunctions,  the 
beam  density  function  can  be  expanded  in  terms  of  the  transverse  eigenfunctions, 
r±),  i.e., 

PM)  =  E  Pmn  (z,t)i/)mn( TjJ  (3.33) 

m,n 

where  pmn(z,t)  are  the  expansion  coefficients,  which  depend  on  the  longitudinal  co¬ 
ordinates  and  time,  m  and  n  are  mode  numbers  of  transverse  eigenfunctions,  and  rj_ 
is  the  transverse  coordinates.  From  Eqn.  (3.33),  there  needs  to  be  enough  transverse 
eigenmodes  in  order  to  resolve  the  transverse  profile  of  the  fields.  The  eigenfunction, 
V;mn(r±),  is  determined  by  the  given  transverse  beam  distribution. 
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For  the  cylindrically  symmetric  beam,  the  azimuthal  mode  number,  m,  is  zero, 
hence,  only  the  radial  eigenmode  number  is  important  for  the  expansion.  For  this 
cylindrical  beam  configuration,  the  eigenfunctions  can  be  simply  expressed  as 


^On(ri) 


1  Mjonr/a ) 

a.y/n  \Mjon)\ 


(3.34) 


where  jon  is  the  nth  root  of  the  zeroth  order  Bessel  function.  In  the  following  sections, 
we  will  use  the  radially  parabolic  shape  of  the  beam  distribution: 

p(r,t)  =  ^0(r6  -  r)  fl  -  S[z  -  z"(t )]  (3.35) 

wrb  \  rb  J 

This  function  has  no  azimuthal  dependence,  i.e.,  it  is  cylindrically  symmetric.  The 
delta  function  indicates  the  beam  has  zero  thickness  in  the  longitudinal  direction. 
z"(t)  contains  the  trajectory  information  of  the  charge,  the  beam  has  total  charge  of 
Q.  Therefore,  the  beam  distribution  function  can  be  expanded  as 


4o(^'r,a) 

yrr  r  *  J  ii 


™b  „  ion  1^1  (ion)  | 2  [jonn/a 


-Jitionn/a)  -  Mjonn/a ) 


S[z-z"(t)]  (3.36) 


The  Eqn.  (3.36)  is  very  similar  to  the  our  Green’s  function  solution  of  the  transverse 
electromagnetic  fields  with  the  given  density  function,  Eqn.  (3.33). 

Figures  3.4  shows  plots  of  the  radial  beam  density  distribution  for  different  eigen¬ 
mode  summations  with  two  different  sizes  of  the  beam  radius,  r^/a  =  0.02433  and 
0.5.  By  comparing  two  Figures  3.4(a)  and  3.4(b),  we  note  that  the  convergence  rates 
of  the  series  of  the  density  function  (Eqn.  (3.36))  strongly  depends  on  the  transverse 
size  of  the  beam.  The  convergence  rate  is  much  more  slow  for  the  small  radial  beam 
size.  For  the  beam  radius  of  r^ja  =  0.5,  even  100  eigenmode  is  enough  to  resolve  the 
density  function.  For  the  beam  radius  of  r^/a  =  0.02433,  the  density  function  is  oscil¬ 
lating  around  the  exact  solution  (solid:red),  therefore,  and  at  least  2,000  eigenmodes 
are  needed  to  get  less  than  1%  numerical  error. 
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(a)  rb/a  =  0.5 


(b)  rb/a  =  0.02433 


(c)  rb/a  =  0.0111 


Figure  3.4:  The  plots  of  the  beam  distribution  function  with  dif¬ 
ferent  beam  sizes.  The  original  function  is  plotted  without  the  ex¬ 
pansion,  and  rmn  is  the  radial  eigenmode  number. 


3.2  Numerical  Implementation 


38 


For  numerical  calculations  of  the  fields,  we  are  required  to  perform  the  transverse 
eigenmode  summations,  a  time  integration,  spatial  integrations  over  three  dimen¬ 
sional  coordinates,  and  the  temporal  and  spatial  derivatives  for  calculating  the  fields. 
In  order  to  accurately  simulate  the  space-charge  fields  within  less  than  1%  error, 
however,  we  need  to  consider  the  computational  requirements  for  these  numerical 
implementations.  Since  the  zero-thickness  parabolic  beam  is  assumed  with  the  cylin- 
drically  symmetric  cavity  geometry  in  this  work,  the  spatial  integrations  are  handled 
analytically.  Due  to  azimuthal  symmetry,  there  is  only  radial  eigenmode  summations 
and  no  azimuthal  mode  summations.  We  now  show  the  computational  requirements 
for  high-accuracy  space-charge  field  calculations. 


3.2.2  Eigenmode  Summation 

With  the  cylindrically  symmetric  beam  distribution  function,  described  in  the  Equa¬ 
tion  (3.35),  the  radial  space-charge  held  is  given  by 

2 


p,  ,n  _  4Q  ^  JoUonr/a) 

r[  ’  _7re0rb2  V  Jonl^Oon)!2 


jonH/ a 


Ji(jonn/a)  -  Mjon‘rb/a ) 


(3.37) 


x  /  dt'  [J0  (jonA_/a)  @(A2 )  -  J0  (j0n A+/a)  0(A2  ) 


where  A^_  are  given  by 


A±2  =  c2  (t-t'f-  [z±z"(t')\ 


"Ul\ l2 


(3.38) 


We  can  estimate  the  required  number  of  radial  modes  to  model  for  a  charge 
density  with  a  form  given  by  Eq  (3.35).  In  general,  a  large  number  of  radial  modes 
are  necessary  to  reconstruct  the  type  of  charge  density  for  a  given  beam  radius,  r^. 
For  large  mode  numbers,  i.e.,  jonr/a  1, 

Jo(jonr/a )  «  J  — —  cos  (j0n-  -  p)  (3.39) 

V  VT j0nr  V  a  4/ 
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is  the  asymptotic  form  for  the  radial  mode  expansion. 

Since  a  cosine  function  has  a  period  2ir,  it  is  necessary  that  in  order  to  have  a 
sufficient  number  of  modes  to  model  the  beam  charge  density, 

joM—  >  27t  (3.40) 

a 

where  M  is  the  largest  required  radial  mode  number.  The  number  of  transverse 
eigenmodes  necessary  for  accurately  determining  the  fields  is  inversely  proportional  to 
the  transverse  size  of  the  beam.  The  necessary  mode  number  is  inversely  proportional 
to  the  ratio  between  the  beam  size,  r b  and  the  radius  of  the  cavity,  a,  i.e., 

27T 

joM  >  — i-  ~  12.57  >  j0, 4,  (3.41) 

n/a 

for  n/a  =  0.5.  Choosing  M  ~  40  satisfies  this  condition.  From  the  inequality  relation, 
The  appropriate  transverse  eigenmode  numbers  can  be  estimated.  For  example,  with 
n/a  =  0.02433,  it  is  necessary  that, 

2n 

Jom  >  —r  ~  258.25  >  jo, 82-  (3.42) 

n/a 

Choosing  M  >  1000  satisfies  this  condition. 

In  Figure  3.5,  the  radial  electric  fields  are  plotted  as  a  function  of  r  for  100, 
500,  1000,  and  2000  eigenmodes  for  the  charge  density  shown  in  Figure  3.4(b)  For 
ct/a  =  0.25  and  z/a  =  0.225,  the  radial  electric  fields  are  oscillating  with  r  for  small 
values  of  eigenmodes.  In  order  to  model  the  fields  within  1%  accuracy,  it  is  necessary 
to  sum  over  at  least  2000  modes. 

3.2.3  Numerical  Time  Integration 

In  order  to  resolve  the  fields  related  with  the  beam  motion,  the  integration  time 
step,  At'  must  be  smaller  than  the  characteristic  time  scale  of  the  system.  Since  the 
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Figure  3.5:  Er  vs.  r  for  different  eigenmode  number,  rmn 
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electron  beam  is  accelerated  rclativistically  by  the  external  rf  field,  the  characteristic 
time  of  the  photoinjector  system,  T^yn  is  given  by 

Tdyn  ~  7 bmece/E.  (3.43) 

where  7 &  is  the  relativistic  Lorentz  factor  of  the  beam.  This  yields  At'  -C  T&yn.  In 
addition  to  this  limit,  however,  we  need  more  restrictions  for  the  time  integration  step 
to  reduce  the  numerical  errors.  There  is  another  restriction  on  the  time-integration. 

The  formula  for  the  space-charge  fields  using  the  Green’s  function  method  con¬ 
tains  time  integral  factors,  which  store  the  trajectory  history  of  the  bunch.  Due  to 
the  complexities  of  the  integrands,  however,  it  is  necessary  to  integrate  these  nu¬ 
merically,  and  to  consider  the  relation  between  the  integration  time  step  and  the 
eigenmode  number.  Each  field  equation  includes  similar  Bessel  function  type  argu¬ 
ments:  Jo(jonA-/a)@(A2_)  and  </oC/onA+/a)0(A+).  The  former  represents  the  real 
charge,  while  the  latter  reads  the  induced  image  charge  by  the  cathode  geometry. 
The  function  \±(z,t]t')  are  given  by 

A|  =  c2(t  -  t’)2  -[z-  z"(t')}2  (3.44) 

where  if  is  the  time  integral  variable  and  varies  from  -00  to  t.  The  step  functions 
represent  the  causality  conditions  of  the  field,  which  occur  due  to  the  finite  speed  of 
light.  For  modeling  the  transient  effects  near  the  front  of  the  bunch,  these  factors 
become  very  important  [21]. 

Figure  3.6  shows  Jo(jonA-/a)0(A?_)  as  functions  of  t1  for  different  transverse  mode 
number,  n,  when  t  is  fixed.  As  n  increases,  Jo(jon\~ / a)Q(\2_)  oscillates  with  time, 
and  these  oscillations  are  smooth.  But  there  are  discontinuities  at  certain  values  of 
t  and  z,  which  are  generated  by  the  step  function.  These  discontinuities  generate 
numerical  time  integration  errors,  so  that  the  time  step,  At',  must  be  small  enough 
to  reduce  errors  around  these  points.  This  imposes  a  maximum  limit  on  the  time 
step. 


3.2  Numerical  Implementation 


42 


(a)  JoOonAl)0(A2_) 


(b)  Jo(io„A+)0(A^) 


Figure  3.6:  Plots  of  the  arguments  of  the  time  integration  for  1 
(red),  10  (green),  100  (blue),  500  (pink),  1000  (light  blue),  and  2000 
(yellow)  eigenmode  numbers. 
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The  time  step  must  be  chosen  sufficiently  small  such  that 


cAt' 

Jo n - <<  2vr, 

a 

so  that  time  integration  over  the  Bessel  functions  are  accurate, 
we  find 


(3.45) 

Using  jo >  2vr, 


At1  « 


0.01rfc 


(3.46) 


From  Figure  3.6,  however,  the  oscillation  time  is  very  short  for  the  large  eigenmode 
number,  thus  we  need  much  smaller  time  step. 


Figure  3.7:  Er  vs.  r  for  different  numerical  integration  time  steps: 
dr'  =  cdt' /a  =  10-4,  10~5,  and  10-6 


Figures  3.7  and  3.8  show  Er  and  Ez  as  functions  of  2  for  different  time  steps,  At', 
with  t/a  =  0.25  Here,  r'  is  the  normalized  time,  ct' /a.  The  maximum  mode  number 
is  chosen  to  be  2000  to  see  only  errors  by  the  numerical  integration  time  step.  For 
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z/a 


Figure  3.8:  Ez  vs.  z  for  different  numerical  integration  time  steps: 
dr'  =  cdt' /a  =  10-4,  10~5,  and  10-6 


3.2  Numerical  Implementation 


45 


the  time  step,  cdt'/a  =  10-4,  10~5andl0-6,  there  is  no  significant  difference  in  Er  vs. 
r.  Since  there  is  no  radial  dependence  in  the  time  integration  factor,  the  effect  of  the 
time  step  size  is  very  small  with  fixed  z  and  t  for  Er. 

For  the  Ez  calculations,  however,  the  longitudinal  field  strongly  depends  on  the 
beam  trajectory,  which  is  the  function  of  the  time  integral  variable,  t' .  Figure  3.8 
shows  the  dependence  of  the  time  step  size  along  the  axis  of  the  symmetry.  The 
oscillation  periods  are  determined  by  the  characteristic  spatial  step  size,  Ax;,  the 
transverse  eigenmode  number,  n,  and  the  integration  step  size,  At'.  The  smaller  Ax; 
reduces  the  period  of  the  oscillation,  however  it  also  requires  much  smaller  At'.  The 
smaller  At'  reduces  the  amplitude  of  the  oscillation.  This  is  generalized  by, 

cAt'  <  Az  -C  rb/'yb  (3.47) 

The  relativistic  Lorentz  factor  of  the  beam,  is  assigned  to  correct  the  relativistic 
effect  of  the  beam.  We  can  see  that  cAt'/a  =  10”4  is  not  small  enough  to  resolve  the 
fields,  even  though  it  is  about  factor  of  200  times  smaller  than  77,/a. 

3.2.4  Multi  Slice  Bunch  Simulation  Method 

Since  the  bunched  beam  has  finite  size  in  the  longitudinal  direction,  the  zero-t hickness 
beam  is  not  adequate  for  designing  the  real  beam.  Thus  the  multi-sliced  bunch  model 
is  introduced  to  make  it  realistic.  We  assume  that  (1)  there  are  a  finite  number  of 
zero-t  hickness  sliced  beams.  (2)  Each  slice  has  an  equal  number  of  particles.  And  (3) 
these  slices  are  uniformly  distributed  over  the  bunch  length.  We  try  to  find  out  how 
many  bunches  are  needed  to  express  the  finite  thickness  beam. 

Suppose  that  the  total  bunch  length  is  ti  =  9  ps.  The  bunch  length  varies  with  the 
frequency,  injection  phase,  and  wave  number  of  the  rf  electric  field.  Since  each  slice  is 
emitted  at  a  different  time,  its  injection  phase  is  slightly  different.  However,  a  9  psec 
bunch  length  is  relatively  short  compared  to  the  time  scale  of  the  rf  field.  Thus,  we 
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can  assume  that  all  slices  have  the  same  trajectory  configurations,  i.e. ,  all  slices  are 
launched  at  the  same  phase,  and  they  are  distributed  uniformly  in  time.  Therefore, 
if  the  central  slice  of  the  bunch  is  at  time  t  =  to,  then  the  ith(l  <  i  <  N/2)  slice  at 
the  head  of  the  bunch  is  at  time  ti  —  to  +  iti/N,  since  the  time  interval  between  two 
adjacent  slice  for  the  N  +  1  slices  is  tjN . 


z/a 


Figure  3.9:  Er  vs.  z  at  r  =  for  different  numbers  of  slices 


Now,  we  set  up  the  uniformly  spaced  (in  time),  equally  charged  N+l  slices,  which 
form  the  complete  bunch.  In  this  section,  we  use  the  AWA  rf  photoinjector  gun  of  the 
ANL.  In  order  to  reduce  other  considerable  errors,  the  following  numerical  parameters 
are  used:  the  eigenmode  number,  n  =  1000,  the  spatial  step  size,  A  z/a  =  10-4,  and 
the  numerical  time  integration  step,  At'  =  5  x  10~5. 

Figures  3.9  to  3.11  show  plots  of  the  space-charge  electric  fields  as  a  function  of 


Ez(47ce0a2/Q)  Ez(47te0a2/Q) 
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(a)  31  slices 


z/a 

(b)  101  slices 


(c)  151  slices 


(d)  301  slices 


Figure  3.10:  Ez  vs.  z  at  r  =  0  for  different  numbers  of  slices 
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(a)  31  slices  (b)  101  slices 


(c)  151  slices 


(d)  301  slices 


Figure  3.11:  Ez  vs.  z  at  r  =  for  the  different  number  of  slices 
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the  normalized  longitudinal  coordinate,  z/a ,  for  the  different  number  of  slices.  The 
fields  are  at  ct/a  =  0.185,  and  those  are  plotted  around  the  bunch. 

In  the  Figure  3.9,  the  radial  held  of  31-slice  model  is  rapidly  oscillating  along  the 
axis,  although  the  errors  are  still  small.  These  oscillations  are  due  to  the  peak  of  each 
slice.  Thus  31  slices  are  not  enough  to  resolve  the  radial  electric  held.  But,  for  more 
than  101  slices,  the  helds  become  smoothed  out  with  small  discrepancies,  which  are 
ignorable  errors. 

In  Figures  3.10  and  3.11,  beyond  the  101  slices,  the  peak  helds  do  not  vary  much 
with  the  number  of  slices.  However,  the  axial  helds  are  oscillating  even  with  151 
slices.  The  amplitudes  of  the  oscillation  are  reduced  with  the  position  of  r.  From 
these  plots,  we  can  estimate  that  more  than  300  slices  can  represent  the  finite  size  of 
the  beam  for  the  bunch  length  of  9  ps.  If  the  bunch  length  is  shortened  or  lengthened, 
the  required  number  of  slices  are  decreased  or  increased,  respectively. 

3.2.5  Benchmark  Modeling  and  Comparison 

It  is  essential  when  developing  a  numerical  algorithm  to  find  benchmarks  by  which 
the  algorithm  may  be  compared.  The  benchmark  that  we  now  present  is  based  on 
an  analytical  result.  In  particular,  we  test  the  case  of  a  source  with  cylindrical  pipe 
of  radius  a  whereby  a  bunch  of  total  charge  Q  is  emitted  from  the  cathode  at  time 
t  —  0  with  a  constant  uniform  velocity  V  =  Vez.  The  bunch  is  assumed  to  be 
cylindrically  symmetric,  have  zero  longitudinal  length,  and  have  finite  size  in  the 
transverse  directions.  We  assume  that  the  radial  density  of  the  bunch  is  parabolic 
with  radius  r>,  and  fixed  in  time.  The  charge  and  current  densities  are  given  by 

p(r,t)  =  p±(r±)S(z  -Vt),  (3.48a) 

and 


Jz  (r,t)  =  p_ l  (rj_)  VS(z-Vt)  , 


(3.48b) 
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where  p±  (rj_)  =  (2Q/nr/Q  {rb  —  r)  (1  —  r2  jr\). 

Now  we  compare  the  fields  of  this  system  to  the  fields  of  another  system  where 
there  are  two  bunches  Q\  —  Q  and  Q2  =  —  Q,  with  velocities  V i  =  Vez  and  V2  = 
— Vez,  respectively.  The  trajectories  of  the  two  bunches  are  configured,  such  that  at 
time  t  =  0,  the  two  bunches  will  overlap.  The  charge  and  current  densities  of  this 
second  system  are  given  by 


pi  (r,  t)  =  (rj.)  S  (z  -  Vt) ,  (3.49a) 

Jzi  (r,t)  =  p±  (rj.)  VS  (z  -  Vt) ,  (3.49b) 

p2(r,t)  = -p±(r±)  8  (z +  Vt) ,  (3.49c) 

and 

JZ2  (r,  t)  =  p_ l  (r_i_)  VS(z  +  Vt).  (3.49d) 

Figure  3.12  shows  the  schematic  comparison  diagrams  between  the  benchmark 
model  and  our  model.  In  the  benchmark  model,  there  is  no  cathode  boundary.  How¬ 
ever,  there  are  two  beams,  which  having  opposite  polarities  and  moving  in  opposite 
direction,  to  represent  the  image  charge  effect  by  the  cathode.  Assuming  that  the 
beam  is  moving  with  the  constant  speed,  e.g.,  v  =  0.9c. 

The  fields  for  each  bunch  in  the  second  system  can  be  easily  computed  by  finding 
the  electric  held  in  the  rest  frame  of  the  bunch,  and  then  Lorentz  transforming  back 
to  the  laboratory  frame.  In  the  second  system,  the  bunch  Q2  represents  the  image 
bunch  which  is  found  in  the  first  system  due  to  the  presence  of  the  cathode  at  z  =  0. 
For  specific  regions  of  space  and  time,  e.g.,  before  wave  reflection  from  the  pipe  has 
occurred  (  t  <  2a/ c  for  r  =  0)  and  after  the  initial  electromagnetic  shock  front  due  to 
bunch  emission  has  passed  ( z  <  ct  for  r  —  0),  the  electromagnetic  fields  produced  by 
the  first  system  will  be  exactly  the  same  as  those  of  the  second  system.  We  note  that, 
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(a)  IRPSS  Model 


-Q 

v 


Q 

v 


(b)  Benchmark  Model 

Figure  3.12:  (a)  IRPSS  simulation  of  a  disk  bunch  of  charge  emit¬ 
ted  at  time  t  —  0  from  the  cathode  surface  moving  uniformly  with 
speed  v  (b)  Analytical  model  of  two  disks  of  charge  moving  uniformly 
in  opposite  directions  for  all  time  and  intersecting  at  t  —  0 
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Figure  3.13:  Plots  of  Er  vs.  r  at  the  bunch  location  z/a  =  0.225, 
showing  the  numerical  scheme  (solid  line)  and  the  benchmark  exam¬ 
ple  (dashed  line). 


in  the  first  system,  the  fields  should  be  exactly  zero  when  z  >  ct,  since  the  causality 
condition  prevents  any  waves  from  propagating  beyond  z  =  ct. 

Figures  3.13  and  3.14  show  plots  of  the  normalized  radial  and  longitudinal  electric 
fields  using  the  beam  densities  given  in  Eqn.  (3.48)  and  the  field  equations  given  by 
Eqns.  (3.18)  and  (3.23).  In  the  simulation,  the  following  parameters  were  chosen:  the 
beam  velocity  is  V  =  0.9c,  the  beam  radius  corresponds  to  the  experimental  beam 
radius  in  the  BNL  2.856  GHz  photocathode  experiment  [40],  i.e.,  r^/a  =  0.0243,  and 
the  time  is  t  —  0.25 a/c.  In  order  to  accurately  compute  the  fields,  we  set  M  =  2000 
and  At'  =  10_6a/c. 

Figure  3.13  shows  a  plot  (solid  line)  of  Er  vs.  r  at  the  bunch’s  location  z/a  =  0.225. 
In  addition,  a  plot  of  the  normalized  radial  electric  field  from  the  second  (benchmark) 
system  is  shown  (dashed  line),  which  uses  the  densities  in  Eqns.  (3.49)  to  compute 
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Figure  3.14:  Plots  of  Ez  vs.  z,  showing  the  numerical  scheme  (solid 
line)  and  the  benchmark  example  (dashed  line). 


the  fields.  From  Figure  3.13,  it  is  seen  that  there  is  excellent  agreement,  i.e. ,  less  that 
1%  error,  between  the  numerical  scheme  and  the  benchmark  example  for  r/a  <  0.085. 
The  difference  between  the  fields  in  Figure  3.13  for  r/a  >  0.085  is  attributed  to  the 
initial  electromagnetic  pulse  generated  at  bunch  emission.  Essentially,  one  can  view 
this  transient  electromagnetic  pulse  as  being  clue  to  transition  radiation,  which  occurs 
when  the  electron  bunch  is  spontaneously  generated  at  the  cathode  surface  at  time 
t  =  0.  Obviously,  this  electromagnetic  pulse  would  not  be  present  in  the  benchmark 
example  since  the  two  bunches  defined  by  Eqns.  (3.49)  are  uniformly  moving  for 
all  time.  It  is  therefore  necessary  to  wait  a  sufficient  amount  of  time  for  the  initial 
transient  pulse  to  pass  in  order  for  the  fields  in  the  benchmark  case  and  the  numerical 
case  to  match. 

Figure  3.14  shows  a  plot  of  Ez  vs.  z  at  r  =  0  for  both  the  numerical  scheme  (solid 
line)  as  well  as  the  benchmark  example  (dashed  line).  We  see  excellent  agreement, 
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i.e.,  less  than  1%  error  between  the  fields  for  locations  which  are  not  close  to  the 
second  pulse  at  z  —  ct.  For  z  >  ct ,  the  numerical  scheme  correctly  shows  that  the 
fields  are  zero  according  to  the  causality  condition. 

3.3  Discussion 

In  the  future,  we  are  planning  to  expand  the  analysis  in  this  section  by  including  the 
effect  of  one  or  more  metallic  irises  on  space-charge  fields.  Irises  are  typically  found 
in  photoinjector  systems,  and  the  first  iris  along  with  the  cathode  and  cavity  walls 
usually  define  the  boundary  of  the  first  half-cell  in  a  photoinjector.  The  irises  will 
affect  the  space-charge  fields  with  the  presence  of  additional  image  charges  and  image 
currents.  There  are  a  variety  of  methods  which  can  be  utilized  for  including  the  effects 
of  irises.  For  example,  in  Ref.  [25,  26],  the  authors  expanded  the  space-charge  fields  in 
each  region  of  injector  space,  i.e.,  before  and  after  the  iris,  using  eigenfunctions  which 
locally  satisfy  the  correct  boundary  conditions.  The  coefficients  of  the  field  expansions 
were  found  numerically  by  applying  boundary  conditions,  such  as  field  continuity,  at 
the  location  of  the  iris.  Usually,  this  technique  results  in  determining  the  elements  of  a 
large  N  xM  matrix  where  N  and  M  are  the  number  of  eigenmodes  used  in  expanding 
the  fields  before  and  after  the  iris.  Another  technique  with  which  the  present  authors 
are  actively  pursuing  is  a  method  developed  by  H.  Bethe  that  expands  the  fields  due 
to  the  iris(es)  using  a  perturbative  multipole  field  expansion  [41].  This  technique 
starts  by  solving  the  space-charge  fields  assuming  that  no  iris  is  present,  i.e.,  that  the 
geometry  of  the  conductor  is  a  pillbox.  These  lowest  order  electromagnetic  fields  can 
then  be  used  to  compute  the  induced  electric  and  magnetic  multipole  moments  of  the 
iris  to  all  orders  in  b/a,  where  the  iris  radius  is  b  and  the  cavity  radius  is  a.  The  total 
space-charge  field  can  then  be  represented  as  the  zeroth  order  fields  plus  the  fields 
due  to  the  electric  and  magnetic  multipoles  for  each  iris. 
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Chapter  4 

Application  to  the  RF 
Photoinjector  Experiments 


In  order  to  simulate  the  beam  dynamics  and  electromagnetic  space-charge  fields  in 
the  rf  photocathode  gun,  including  those  near  the  cathode,  a  new  code  called,  IRPSS 
(Indiana  Rf  Photocathode  Source  Simulator),  has  been  developed.  IRPSS  can  cal¬ 
culate  the  electromagnetic  space-charge  fields  based  on  a  time-dependent  Green’s 
function  method  [42,  21,  43,  44],  and  it  can  track  the  beam  evolution  throughout  the 
first  half  cell.  In  order  for  a  one-to-one  comparison,  we  have  written  an  electrostatic 
code,  which  is  also  based  on  a  Green’s  function  method.  To  be  specific,  the  elec¬ 
trostatic  code  incorporates  both  E  and  B  space-charge  fields,  but  does  not  include 
time-retardation  effects.  While  each  code  has  a  different  space-charge  solver,  both 
codes  utilize  identical  tracking  algorithms. 
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4.1  Simulation  Set-Up  for  the  RF  Photoinjector 


4.1.1  Prescribed  Trajectory 


In  the  IRPSS  code,  we  first  calculate  the  electromagnetic  space-charge  potentials 
and  fields  for  a  given  charge  and  current  densities,  which  satisfy  the  continuity  equa¬ 
tion.  The  space-charge  fields  are  computed  using  the  time-dependent  Green’s  function 
method.  With  these  self-fields,  we  track  the  motions  of  the  bunched  beam  in  the  rf 
photoinjector  under  the  influence  of  external  fields,  such  as  rf  and  solenoid  fields. 

IRPSS  constructs  the  charge  and  current  densities  of  the  beam  using  slices,  i.e., 
zero  thickness  disks.  In  particular,  the  choice  of  p  (r)  and  Jz  (r)  is  of  the  form 

N 

P(r)  =  ^2ai(r)SlZ  ~  zi(t)\  >  (4'la) 

i=  1 


and 


N 


Jz(r)  =J2ai(r)S[z-zi(t)\ 


i— 1 


dz"(t) 

dt 


(4.1b) 


where  Uj(r)  is  the  charge  per  unit  area  of  the  ith  slice,  z"{t)  is  the  longitudinal  location 
of  the  ith  slice,  and  N  is  the  number  of  slices  on  the  simulation.  The  required  number 
of  slices  to  form  the  bunch  was  described  in  the  previous  Chapter  3. 

In  order  to  illustrate  this  numerically,  we  use  the  ANL  AWA  1.3  GHz  gun  parame¬ 
ters  as  listed  in  Table.  3.1.  We  make  the  simplifying  assumption  that  the  trajectories, 
z"(t),  are  specified  by  the  external  rf  field.  Physically,  one  could  view  this  simulation 
as  the  space-charge  fields  being  small  compared  to  the  rf  fields.  The  trajectory  of 
each  of  the  slices  can  be  computed  from  the  relativistic  equation  of  motion  for  the 
slice  in  the  rf  field  [11]. 


dzT_ 

dt 


(4.2a) 
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ct/a 


Figure  4.1:  The  trajectory,  velocity,  and  acceleration  of  a  particle 
in  ANL  AWA  1.3  GHz  rf  photoinjector  are  plotted  as  functions  of 
r  =  ct/a  with  the  rf  peak  field  of  77  MV/m. 
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and 


dP * 
dt 


eE?{r,z,t), 


(4.2b) 


where  me  is  the  mass  of  an  electron  and  ErJ(r,  z,  t )  is  the  external  rf  field.  Figure  4.1 
shows  the  trajectory,  velocity,  and  Lorentz  factor,  7,  of  a  particle  in  the  ANL  AWA 
1.3  GHz  rf  photoinjector.  To  calculate  the  trajectories,  the  TM0i  mode  rf  field  is 
used. 


4.1.2  External  Fields 

In  the  rf  photoinjector,  the  electron  beam  emitted  from  the  cathode  is  accelerating 
by  the  external  rf  held.  The  rf  held  accelerates  the  beam  from  rest  to  a  relativistic 
velocity.  There  are  numerous  codes  for  modeling  the  rf  cavities  and  for  computing  the 
holds  in  those  systems,  such  as  SUPERFISH  [33],  or  MWS  (Micro  Wave  Studio)  [47]. 
In  our  simulation,  we  calculate  the  rf  helds  using  SUPERFISH,  and  the  held  data  is 
easily  imported  into  our  simulation  code,  IRPSS.  Figure  4.2  shows  the  external  rf  and 
solenoid  helds  which  have  been  used  in  our  beam  dynamics  simulation.  Those  helds 
are  computed  using  parameters  of  the  ANL  AWA  1.3  GHz  rf  photoinjector,  listed  in 
Table  3.1.  In  the  following  numerical  examples  for  the  beam  dynamics,  we  will  use 
the  rf  held  maps  calculated  using  SUPERFISH. 

In  SUPERFISH,  the  rf  system  is  assumed  to  be  cylindrically  symmetric.  Fully 
3-dimensional  rf  held  maps  with  complicated  geometries  can  be  calculated  with  other 
simulation  codes,  e.g.,  MWS.  In  IRPSS,  the  external  held  maps  are  imported  in  3- 
dimensional  coordinates,  hence  we  can  use  not  only  rf  held  maps  of  SUPERFISH,  but 
also  those  of  MWS.  Since  our  goal  is  to  see  the  zeroth  order  symmetric  space-charge 
effect,  the  geometry  of  the  rf  gun  is  assumed  to  be  the  axially  symmetric,  and  this  is 
a  good  approximation. 

Alternatively,  one  can  use  TM0io  mode  rf  helds,  which  is  mostly  dominant  through 
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Figure  4.2:  The  external  rf  (dark  blue)  and  solenoid  (green)  fields 
are  plotted  in  the  schematic  diagram  (light  blue)  of  the  ANL  AWA 
1.3  GHz  rf  photoinjector.  For  comparison,  the  rf  fields  (orange)  cal¬ 
culated  from  SUPERFISH  and  those  of  TM0io  mode  are  presented. 
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the  entire  rf  gun.  The  advantage  of  TMoio  mode  rf  fields  is  that  the  analytical  forms 
are  well  known  and  ready  to  use.  The  TM0io  mode  electromagnetic  rf  fields  in  the 
pillbox  cavity  with  sin  convention  are  given  by  [6,  39] 


Ez 

Er 

Be 


—E0  cos (kzz)  sin(c jt  +  <po), 
E0rkz 


2 

E0rkz 

2c 


sin^z)  sin(o;t  +  (fo ), 
cos (kzz)  cos(u it  +  (po), 


(4.3) 


where  Eq  is  the  peak  rf  field,  kz  is  the  longitudinal  wave  number,  c o  is  the  angular  rf 
frequency,  and  (po  is  the  injection  rf  phase.  The  vector  and  scalar  potential  can  be 
found  from  the  following  three  relations: 

(1)  E  =  -V»  -  §  , 

(2)  B  =  V  x  A  , 

(3)  V  •  A  +  =  0  (Lorentz  gauge  condition). 

Then,  the  potentials  are  given  by 


$  : 

Ar 

Az 


_Eo_ 
2  kz 

E0r 

2c 

Eo 

2  kzc 


sin  (kzz) 

i  sin  (cut  +  <p0), 

(4.4a) 

sin  (kzz) 

cos  (cut  +  <po), 

(4.4b) 

cos [kzz] 

1  cos(c dt  +  (po). 

(4.4c) 

The  solenoid  magnetic  fields  in  IRPSS  simulation  are  also  imported  from  the 
external  field  maps.  The  magnetic  fields  can  be  computed  using  POISSON  [33]  or 
MAFIA  [19].  Since  our  simulation  is  concentrated  on  the  first  half  cell  in  the  rf  gun, 
and  the  solenoid  magnetic  field  is  relatively  small  as  <  1  kG,  we  do  not  show  its 
details. 
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4.1.3  Equations  of  Motion 


After  the  full  calculations  of  electromagnetic  space-charge  fields  in  the  first  part  of 
IRPSS,  the  second  part  in  IRPSS  code  tracks  the  motion  of  the  electron  beam  with 
space-charge  and  external  fields.  One  may  simulate  the  beam  dynamics  using  either 
the  Hamiltonian  of  the  system  or  the  Lorentz  force  law.  The  Hamiltonian  of  an 
electron  in  the  cylindrical  coordinates  is  given  by 


H  =  ecf>  +  c 


tti2c2  +  {Pr  —  eAr)  +  I - J  -\-  (Pz  —  gAz) 


1  1/2 


(4.5) 


In  IRPSS,  we  prefer  to  use  the  Lorentz  force  law, 


F  =  q  [E  +  (/3  x  cB)] ,  (4.6) 

since  this  method  does  not  require  additional  derivatives  of  potential  maps.  In  the 
cylindrically  symmetric  system,  we  have  only  the  following  field  components. 

(1)  Space-Charge  fields:  EfG,  Ezc,  and  BgG , 

(2)  Rf  fields:  EFF,  EfF,  and  BfF, 

(3)  Solenoid  fields:  Bfc  and  Rfc. 

The  2nd  order  differential  equations  in  the  Lorentz  force  law  with  these  fields 
are  numerically  simulated  using  the  4th  order  Runge-Kutta  method  [36],  Runge- 
Kutta  methods  have  the  high-order  local  truncation  error  of  the  Taylor  methods  while 
eliminating  the  need  to  compute  and  evaluate  the  derivatives  of  /  (r ,t).  To  achieve 
highly  accurate  results,  CERNlib  is  implemented  using  the  Runge-Kutta  method  [48]. 
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4.2  Space-Charge  Fields  and  Beam  Dynamics 

4.2.1  Calculations  of  Space-Charge  Fields 

In  most  of  rf  simulation  codes,  such  as,  PARMELA,  the  space-charge  fields  are  cal¬ 
culated  electrostatically.  However,  IRPSS  uses  a  fully  electromagnetic  method  for 
the  space-charge  held  calculations.  The  differences  in  held  calculations  can  cause 
differences  in  the  beam  dynamics.  In  this  section,  we  hrst  show  the  physics  of  space- 
charge  helds  near  the  cathode,  and  compare  the  electromagnetic  and  electrostatic 
space-charge  helds. 

In  order  to  compare  the  space-charge  helds  of  IRPSS  with  electrostatic  codes,  we 
simplify  the  photoinjector  system  by  assuming  that  it  has  no  iris  and  the  cavity  is 
cylindrically  symmetric.  Since  the  iris  does  not  have  much  of  an  effect  on  the  beam 
dynamics  for  early  times,  this  assumption  is  an  excellent  approximation,  while  we 
investigate  the  physics  near  the  cathode. 

The  rf  held  for  the  beam  accelerating  in  the  gun  is  calculated  using  SUPER¬ 
FISH  [33],  and  the  peak  rf  Ez  held  is  50  MV/m  with  the  injection  phase  of  65 
degrees.  The  azimuthal  rf  magnetic  held  Bg  is  also  included  in  the  simulation.  Here, 
only  rf  helds  are  used  for  the  external  held  in  order  to  focus  on  the  effect  of  the  space- 
charge  helds.  The  external  solenoidal  helds,  which  are  normally  used  for  omittance 
compensation  [12],  are  very  small  near  the  cathode  and  are  excluded. 

We  use  a  radial  parabolic  shape  for  the  transverse  beam  distribution  and  a  uniform 
longitudinal  distribution  for  these  simulations: 

p(r,z,t)  =5^“W0(u-'r)  5[z~zi(t)\  (4-7) 

i  nr2  V  rbJ 

where  Q  is  the  bunch  charge,  r&  is  the  beam  radius  (1  mm).  For  the  uniform  longitu¬ 
dinal  beam  distribution,  a  multi-slice  method  is  used.  In  particular,  a  bunched  beam 
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is  broken  up  into  a  number  of  zero-thickness  slices.  In  Eq.  (4.7),  the  z"(t)  indicates 
the  location  of  the  ith  slice  with  respect  to  time.  We  assume  that  each  beam  slice 
has  an  equal  amount  of  charge.  To  reduce  the  computational  time  and  errors  and 
to  achieve  accurate  calculations  of  space-charge  fields,  an  adequate  number  of  zero 
thickness  slices  are  uniformly  distributed  within  the  desired  bunch  length.  In  the 
space-charge  held  calculation,  we  choose  a  bunch  length  of  1.2  ps.  For  this  bunch 
length,  41  longitudinal  slices  are  enough  to  smoothly  model  the  space-charge  fields. 
In  Chapter  3,  we  showed  the  computational  requirements  for  IRPSS,  such  as  the 
minimum  eigenmode  numbers  and  the  minimum  numerical  integration  time  step  for 
achieving  highly  accurate  space-charge  held  calculations. 

In  IRPSS,  the  electromagnetic  space-charge  fields  are  solved  from  Maxwell’s  equa¬ 
tions  for  conductor  boundary  conditions.  These  solutions  can  be  obtained  using  a 
time-dependent  Green’s  function  method.  For  cylindrically  symmetric  pipe  bound¬ 
ary  conditions  with  a  hat  cathode,  all  helds  can  be  computed  analytically  and  are 
given  by, 


Er(r,z,t)  =  ^ 

n= 0 
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and 
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(4.8c) 


where  ft,  —  ri,/a,f  —  r/a,  z  —  z/a,  r  —  ct/a,  t'  —  ct' /a,  and  A2  =  (r  —  r')2  — (i  ±  z')2 
For  the  cylindrical  cavity,  ipmn  are  the  transverse  eigenfunctions.  In  general,  spa¬ 
tial  and  time  integrations  can  be  performed  numerically.  However,  in  this  simulation, 
we  specify  p  in  Eqn.  (4.7)  for  the  transverse  direction,  and  use  zero  thickness  slices 
longitudinally  to  model  the  beam. 

In  addition  to  IRPSS,  we  have  also  written  an  electrostatic  code  using  a  time- 
independent  Green’s  function  approach.  The  electrostatic  space-charge  potential  of 
the  real  disk  beam  in  the  rest  frame  for  a  cylindrically  symmetric  waveguide  with  a 
flat  cathode  can  be  calculated  as  [34]: 

<Kr)  =  —  /  Gy  (r;  r')  p  (r')  dr',  (4.9) 

eo  J 


where 
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mn  (r±)^nir'±)„  _ 3  On  |  g _ Zq  | 
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’  ’  2  ^ 
mn 


v.Lmn 


(4.10) 


and  a  is  the  cavity  radius  (9.08  cm).  In  order  to  include  the  image  charge  effect, 
the  scalar  potential  can  be  described  as  the  superposition  of  potentials  of  real  and 
image  disks.  Since  the  beam  is  at  rest  in  this  frame,  the  vector  potential,  Az  is  zero. 
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The  scalar  and  vector  potentials  in  the  laboratory  frame  can  be  obtained  by  Lorentz 
transformations,  and  these  yield  the  electrostatic  space-charge  fields: 
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e  o  u  1  —  e 

where  zo  is  the  location  of  a  disk  beam.  The  detailed  derivation  of  Eqn  (4.11)  is  in 
Appendix  B. 

In  order  to  check  the  validity  of  the  electrostatic  code,  we  also  formulated  an 
exact  bunch  charge  and  current  density  fluid  model  assuming  zero  current  density 
compressibility,  which  represents  the  physics  at  beam  emission.  We  have  found  that 
both  the  electrostatic  fluid  model  and  the  electrostatic  code  agree  precisely.  One  of 
the  advantages  of  developing  the  analytical  fluid  model  is  that  we  can  produce  an 
exact  formula  for  the  longitudinal  space-charge  field  at  the  critical  location  r  =  0.0 
m  and  z  =  0.0  m.  This  formula  is  given  by, 


4 Q  ^  / - n - Erf  (VjomZb/a) 

7re°rb  m  V  JOmZb/a  jlmJl(jom) 


(4.12) 
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where  Erf  is  the  error  function,  and  rb  and  zb  are  the  bunch  radius  and  bunch  length, 
respectively.  The  detailed  derivation  of  Eqn  (4.12)  is  in  Appendix  B.2.  When  we 


4.2  Space-Charge  Fields  and  Beam  Dynamics 


66 


compare  the  two  longitudinal  electrostatic  fields  at  r  =  0  and  z  —  0,  which  is  the 
center  of  the  cathode,  as  the  full  bunch  is  emitted,  they  agree  extremely  well. 

Figure  4.3  shows  the  critical  Ez  space-charge  fields  at  r  =  0.0  and  z  =  0.0  for 
different  laser  pulse  lengths  at  0.1  nC  (for  ry  =  1  mm)  and  10  nC  (for  ry  =  1  cm).  The 
primary  result  is  that  as  the  laser  pulse  length  is  increased,  the  discrepancy  between 
electrostatic  and  electromagnetic  space-charge  fields  is  gradually  increased.  The  size 
of  the  discrepancy  also  depends  on  the  beam  radius.  Qualitatively,  as  the  beam 
radius  is  increased,  i.e.,  the  beam  becomes  more  pancake  like,  and  the  discrepancy 
becomes  smaller.  Moreover,  the  maximum  electric  field  for  both  the  electrostatic,  as 
well  as  the  electromagnetic  cases,  occurs  when  they  have  a  zero  bunch  length  beam, 
which  acts  as  a  sheet  of  charge  near  the  cathode.  In  this  case,  the  electric  field  is 
exactly  found  to  be  cr/e0  where  a  is  the  local  charge  density  (charge/area)  on  the 
beam  axis.  For  a  parabolic  (radial)  and  uniform  (longitudinal)  beam  with  Q  =  100 
pC,  a  is  given  by  2Q/neorrl,  hence  a/e o  =  —7.19  MV/m.  Both  electromagnetic  and 
electrostatic  results  show  that  they  agree  extremely  well  for  shorter  bunch  lengths, 
i.e.,  pancake-like  beam,  as  seen  in  Figure  4.3(b). 

Figure  4.4  shows  the  percentage  difference  between  electromagnetic  and  electro¬ 
static  critical  Ez  fields.  And  this  effect  can  become  larger  by  increasing  the  bunch 
length  with  fixed  beam  radius.  However,  the  percentage  discrepancy  between  elec¬ 
tromagnetic  and  electrostatic  fields  has  a  maximum  around  2.4%  for  ry  =  1  mm  with 
a  laser  pulse  length  of  10  psec  and  around  1.6%  discrepancy  for  ry  =  1  cm. 

In  Figures  4.5  and  4.6,  we  compare  electromagnetic  and  electrostatic  space-charge 
fields  for  a  1  nC  bunched  beam.  The  space-charge  fields,  F /Q,  are  computed  from 
the  Lorentz  force  law,  F  =  Q  [E  +  f3  x  (cB)],  where  E  and  B  are  proportional  to  the 
bunch  charge  Q.  Figure  4.5  shows  Fr/Q  vs.  r,  for  various  times.  At  t  =  0.005a/c, 
a  bunch  is  just  launched,  and  there  are  significant  discrepancies  between  the  two 
space-charge  fields.  For  a  large  amount  of  beam  charge,  such  as  1  nC ,  the  beam  feels 
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Laser  Pulse  Length:  tb  (ps) 

(a)  rb  =  1  mm 


> 


Laser  Pulse  Length:  tb  (ps) 

(b)  r;,  =  1  cm 


Figure  4.3:  The  longitudinal  space-charge  fields  at  r  =  0.0  and 
z  =  0.0  computed  by  IRPSS  (red)  and  electrostatic  methods(blue) 
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Figure  4.4:  Percentage  difference  between  critical  electromagnetic 
and  electrostatic  Ez  space-charge  fields  for  different  beam  radii 
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(a)  t  =  0.005a/c 


(b)  t  =  0.020a/c 


(c)  t  =  0.040a/c 


Figure  4.5:  Electromagnetic  (IRPSS)  and  electrostatic  space- 
charge  fields  for  different  times:  Fr/Q  vs.  r  at  the  center  of  the 
beam 
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more  of  a  defocusing  force  with  the  electrostatic  space-charge  fields  in  the  transverse 
direction.  Hence,  this  effect  causes  a  growth  in  beam  phase  space.  These  discrepancies 
in  space-charge  fields  between  the  two  algorithms  become  smaller  with  time.  For  a 
small  amount  of  beam  charge,  it  becomes  more  difficult  to  ford  the  difference  between 
two  fields,  since  the  transverse  space-charge  fields  are  much  smaller  compared  to  the 
transverse  rf  fields. 

Comparisons  of  the  longitudinal  space-charge  fields  are  shown  in  Figure  4.6.  As 
can  be  seen,  the  electromagnetic  Ez  space-charge  held  has  a  more  negative  value  than 
the  electrostatic  Ez  space-charge  held  near  the  cathode  right  after  beam  emission. 
However,  the  discrepancy,  between  electrostatic  and  electromagnetic  space-charge 
holds  in  the  longitudinal  direction,  for  our  example,  is  less  pronounced  than  in  the 
radial  direction,  because  of  the  choice  of  short  bunch  length,  i.e.,  1.2  psec.  The 
longitudinal  discrepancy  becomes  more  apparent  for  longer  bunch  lengths  such  as  in 
the  study  by  A.  Candel  et  al  which  looked  at  a  10  and  15  psec  bunch  length  cases  [23]. 
Like  the  transverse  space-charge  helds,  however,  the  longitudinal  space-charge  helds 
agree  for  higher  time. 

We  could  roughly  estimate  these  discrepancies  by  assuming  the  beam  as  a  thin 
disk.  Since  the  initial  bunch  length  of  1.2  psec  is  relatively  smaller  than  the  beam 
radius  of  ry  =  1  mm,  the  thin  disk  assumption  is  valid.  For  sufficiently  small  times, 
i.e.,  t  «  Tb/c,  the  edge  effects  of  the  beam  for  the  electromagnetic  space-charge  helds 
at  r  —  0  and  z  =  0  is  zero  due  to  the  finite  speed  of  light.  Hence,  the  longitudinal 
space-charge  helds  near  the  cathode  at  r  —  0  resemble  those  of  an  infinite  sheet  of 
charge  at  z"{t)  with  charge  density  cq  and  an  infinite  sheet  of  image  charge  at  —z"{t) 
with  charge  density  —  cr0.  In  this  case,  cr0  =  2Q/'K,rl  which  is  the  charge  density 
at  r  =  0  for  a  parabolic  radial  beam  distribution.  Therefore,  for  short  times,  the 
electromagnetic  longitudinal  held  is  (x0/e0.  In  the  electrostatic  case,  however,  the 
effects  of  the  beam  edges  are  non-negligible  and  give  a  reduction  to  Ez  as  shown  in 
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Figure  4.6:  The  longitudinal  electromagnetic  (IRPSS)  and  electro¬ 
static  space-charge  fields  for  different  times:  Fz/Q  vs.  z  at  r  =  0.0 
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Figure  4.6,  near  the  head  and  tail  of  the  bunch. 

4.2.2  Beam  Dynamics  with  Space-Charge  Effects 

In  the  previous  section,  we  showed  the  calculations  of  the  time-dependent  and  time- 
independent  space-charge  fields  (forces),  and  explained  the  differences  between  IRPSS 
and  the  electrostatic  calculations.  Now,  we  examine  the  beam  dynamics  near  the 
cathode,  where  the  space-charge  effects  are  extremely  important.  As  we  have  seen 
in  the  previous  section,  the  electromagnetic  and  electrostatic  space-charge  fields  have 
disagreements  near  the  cathode,  and  these  are  getting  smaller  with  time.  In  Fig¬ 
ures  4. 7-4. 9,  different  locations  (z  =  0.0  cm  (cathode),  z  =  0.003  cm,  and  z  =  0.03 
cm),  all  within  the  first  cell  and  near  the  cathode,  are  chosen  to  see  differences  of  the 
transverse  phasespaces  ( pr/mc  vs.  r / a)  between  two  space-charge  held  calculations. 

In  order  to  check  the  validity  of  our  tracking  algorithm,  we  first  run  the  simulations 
only  with  rf  holds  (zero  space-charge  helds).  We  compare  these  beam  dynamics  results 
with  those  of  PARMELA,  which  are  also  run  without  the  space-charge  helds.  At  the 
cathode,  the  beam  starts  with  zero  momentum,  i.e.,  zero  thermal  emittance.  Hence, 
particles  are  on  the  r  axis  on  the  transverse  phasespace  plots.  With  only  rf  helds, 
particles  are  radially  focused  near  the  cathode.  Figure  4.7  shows  that  the  two  codes 
agree  very  well. 

For  sufficiently  large  amounts  of  charge,  the  evolutions  of  the  phasespace  are  very 
different.  In  Figures  4.8  and  4.9,  we  include  the  effects  of  space-charge.  Since  the 
initial  beam  distribution  in  the  radial  direction  is  parabolic,  the  radial  space-charge 
forces  are  maximum  at  slightly  less  than  r  =  rb,  as  seen  in  Fig  4.5.  Therefore,  the 
radial  beam  edge  feels  less  defocusing  force  than  the  particles  with  the  maximum 
space-charge  force.  In  phasespace  plots,  we  could  see  these  effects  in  both  IRPSS  and 
the  Green’s  function  based  electrostatic  code. 
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Figure  4.7:  Phasespace  plots  with  zero  space-charge  field 


4.2  Space-Charge  Fields  and  Beam  Dynamics 


74 


r  (mm) 

(a)  z  =  0.003  cm 


Figure  4.8:  Phasespace  plots  for  Q  =  100  pC 
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Figure  4.9:  Phasespace  plots  for  Q  =  500  pC 
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As  we  increase  the  bunch  charge,  the  differences  between  transverse  phasespace 
plots  become  apparent  at  1  pC,  and  then  become  more  significant  for  larger  bunch 
charges.  Since  the  electrostatic  space-charge  fields  near  the  cathode  are  larger  than 
the  electromagnetic  space-charge  fields,  we  see  less  defocusing  in  IRPSS  phasespace 
plots.  As  seen  in  Figures  4.8(a),  and  4.9(a),  IRPSS  and  the  electrostatic  phases¬ 
pace  plots  have  large  discrepancies,  since  particles  with  electrostatic  fields  feel  much 
stronger  defocusing  forces.  In  general,  we  find  a  larger  phasespace  size  in  IRPSS 
simulations  near  the  cathode  compared  to  electrostatic  simulations  due  to  the  larger 
variation  of  Fr  in  IRPSS  with  respect  to  time.  Since  the  two  transverse  space-charge 
fields  are  getting  closer  with  time  as  shown  in  Figure  4.5,  phasespace  plots  for  both 
electromagnetic  (IRPSS)  and  electrostatic  space-charge  fields  look  similar  at  later 
times  as  shown  in  Figures  4.8(b),  and  4.9(b). 

For  a  highly  charged  beam,  the  longitudinal  space-charge  forces  from  the  head 
of  the  bunch  cause  the  tail  of  the  bunch  to  hit  back  onto  the  cathode  causing  beam 
loss.  There  is  no  beam  loss  for  lower  beam  charge,  e.g.  less  than  100  pC,  since  the 
accelerating  Ez  fields  are  dominant.  Beam  loss  starts  to  occur  at  640  pC,  where  the 
space-charge  fields  match  the  rf  fields.  From  Figures  4.6,  we  can  estimate  that  beam 
loss  will  occur  at  a  slightly  smaller  beam  charge  value  for  the  electromagnetic  case 
since  the  electromagnetic  Ez  space-charge  field  at  a  bunch  tail  is  slightly  larger  in 
magnitude  than  the  electrostatic  Ez  space-charge  field.  For  example,  for  1  nC  beam 
at  ry  =  1  mm  and  a  bunch  length  of  1.2  psec ,  IRPSS  expects  about  32.6%  beam  loss, 
while  the  electrostatic  code  yields  32.4%  beam  loss.  These  small  discrepancies  are 
also  found  in  longitudinal  energy  spread  comparisons  which  show  differences  of  less 
than  1%. 

Figure  4.10  shows  comparisons  of  transverse  omittances  between  IRPSS  and  elec¬ 
trostatic  calculations  with  different  bunch  charges.  For  lower  beam  charge,  such  as 
1  pC,  rf  forces  are  dominant  and  space-charge  forces  are  relatively  small,  so  that 
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Figure  4.10:  Normalized  rms  transverse  emittances 
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omittances  are  small  as  shown  in  Figure  4.10(a).  For  Q  =  100  pC ,  however,  the  beam 
feels  much  stronger  space-charge  forces,  hence,  this  causes  larger  omittance  than 
lower  beam  charge  as  shown  in  Figure  4.10(b).  The  discrepancies  in  space-charge 
forces  leads  to  a  difference  in  omittances,  which  is  about  8  ~  9%  at  z  —  0.05  cm. 

4.3  Experimental  Measurements 

We  performed  an  experimental  beam  loss  measurement  on  the  1.3  GHz  rf  gun  at  the 
ANL  AWA  experiment  [45] .  The  main  purpose  of  the  Argonne  Wakefield  Accelerator 
is  for  high  current  beam  generation  and  diagnostics,  and  high-gradient  wakehcld  ex¬ 
periments  [27,  49].  Figure  4.11(a)  shows  the  AWA  rf  photoinjector  and  the  schematic 
diagram  of  its  beamline.  It  consists  of  a  1.5  cell  rf  cavity,  CS^Te  semiconductor  pho¬ 
tocathode,  YAG  laser,  solenoid  magnets,  and  beam  diagnostic  systems.  The  gun  is 
operated  with  50  MV jm  rf  peak  held,  65  degree  injection  phase. 

Before  guiding  the  laser  light  to  the  cathode,  we  measure  the  laser  spot  size  and 
intensity  using  UV  (Ultra-Violet)  camera.  Although  the  efficiency  of  the  photocath¬ 
ode  is  not  100  %  in  generating  the  electron  beam,  we  assume  that  the  laser  intensity 
is  equivalent  to  the  initial  beam  charge.  Experimentally,  it  is  difficult  to  measure  the 
beam  charge  at  the  inside  of  the  gun.  Hence,  we  measure  the  emitted  beam  charge 
by  taking  the  beam  intensity  and  its  profile  through  the  2nd  diagnostic  window  after 
the  gun.  In  order  to  see  the  effect  of  the  bunch  length  in  the  beam  loss  measure¬ 
ment,  we  choose  two  laser  pulse  lengths,  3.4  and  10.4  psec  in  FWHM  (Full-Width 
Half-Maximum).  The  laser  pulse  length  is  equivalent  to  the  bunch  length.  The  laser 
spot  size  is  2  mm.  i.e.,  the  beam  radius  is  1  mm. 

Figure  4.12  shows  plots  of  the  measured  bunch  charge  from  the  exit  of  the  gun 
as  a  function  of  laser  pulse  intensity.  If  no  beam  loss  were  to  occur,  then  the  plot 
should  be  linear  with  a  uniform  slope.  However,  at  a  critical  bunch  charge,  i.e., 
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(a)  ANL  AWA  1.3  GHz  rf  photoinjector 


Cathode 


Beam  dump 


(b)  Schematic  diagram  of  Argonne  Wakefield  Accelerator  beamline  [49] 


Figure  4.11:  Experimental  set  up  for  beam  loss  measurements 


4.3  Experimental  Measurements 


80 


Figure  4.12:  Experimental  results  of  bunch  charge  vs.  laser  pulse 
intensity  for  different  bunch  lengths. 
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E ^  =  £;™ticai ^  fQr  gxec|  ]aser  pulse  length  and  radius,  one  would  expect  beam  loss 
to  occur  and  a  reduction  in  the  slope  of  the  curve.  The  visual  inspection  of  the 
preliminary  data  suggest  changes  in  slopes,  i.e.,  onsets  of  beam  losses,  at  around 
300  —  400  pC  for  3.4  psec  and  500  —  600  pC  for  10.4  psec.  However,  these  are  lower 
than  onr  theoretical  predictions  of  1.2  nC  for  a  uniform  radial  bunch  distribution. 
Possible  explanations  for  these  are  non-uniformity  of  the  laser  spot  size,  effects  due 
to  the  acceptance  in  the  first-half  cell,  and  uncertainties  in  measurements  of  laser 
intensities  and  bunch  charges.  Further  analysis  is  needed  to  understand  differences 
between  prediction  and  experiment.  We  plan  to  continue  these  measurements,  and 
look  at  possible  ways  to  increase  the  Ez  field  discrepancy  by  changing  the  magnitude 
of  the  rf  held  and  injection  phase. 

4.4  Discussion 

IRPSS  is  a  novel  space-charge  solver  and  simulation  code  for  modeling  photoinjectors. 
Since  it  uses  a  Green’s  function  approach  to  calculate  electromagnetic  space-charge 
fields,  it  is  capable  of  simulating  beams  with  arbitrarily  small  bunch  lengths.  Be¬ 
cause  IRPSS  self-consistently  models  the  space-charge  fields  electromagnetically,  ef¬ 
fects  such  as  a  finite  time-retardation  can  be  included.  In  general,  electrostatic  codes, 
such  as  PARMELA  may  overestimate  the  transverse  space-charge  holds  and  underes¬ 
timate  the  longitudinal  space-charge  fields  compared  to  IRPSS  leading  to  differences 
in  beam  phasespace  near  the  cathode  and  beam  loss. 

We  present  the  space-charge  effects  in  the  ANL  AWA  rf  photocathode  gun  near 
the  cathode  by  comparing  the  phasespace  and  omittance  between  IRPSS  and  an 
electrostatic  Green’s  function  code.  In  general,  we  have  found  that  the  phasespace 
plots  and  emittances  predicted  by  IRPSS  and  the  electrostatic  code  can  be  different 
due  to  the  space-charge  fields. 
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Chapter  5 

Transverse  Current  Modeling 


We  generalize  the  exact  formalism  for  the  space-charge  fields  of  a  cylindrically  sym¬ 
metric  beam  in  a  conducting  circular  pipe.  This  formalism  also  includes  the  effect  of 
the  transverse  currents,  which  was  ignored  for  modeling  the  rf  photoinjector  in  the  pre¬ 
vious  chapters.  With  the  series  expansions  of  the  charge  and  current  densities  in  terms 
of  the  Bessel  functions,  the  electromagnetic  space-charge  fields  are  constructed  using 
the  time-dependent  Green’s  function  method  in  the  cylindrical  conducting  boundary 
system.  This  formalism  can  be  used  to  model  the  high  space-charge  systems,  such  as 
high-power  microwave  sources.  Numerical  studies  will  be  shown  as  an  example  of  the 
space-charge  held  calculations  for  a  radially  oscillating  beam. 

5.1  Electromagnetic  Field  for  Circularly  Symmet¬ 
ric  Sources 

In  this  section,  we  derive  the  exact  electromagnetic  space-charge  fields  for  circularly 
symmetric  charged  sources  using  a  time-dependent  Green’s  function  method  [50].  We 
assume  that  the  system  has  a  cylindrically  symmetric  conducting  pipe  with  radius  of 
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a.  The  charge  and  current  densities  are  denoted  by  p(r,  z,  t )  and  J(r,  z,  t)  in  cylindrical 
coordinates.  These  beam  densities  are  only  restricted  by  the  self-consistent  continuity 
equation,  i.e.,  dp/dt  +  V  •  J  =  0. 

5.1.1  Expansions  of  Charge  and  Current  Densities 

First,  we  start  by  making  the  judicious  representations  of  the  charge  and  current 
densities  into  series  of  Bessel  functions.  Since  the  beam  source  and  the  system  are 
cylindrically  symmetric,  Bessel  functions  are  helpful  to  find  the  solutions  of  wave 
equations,  which  define  the  electric  and  magnetic  field  components.  By  noting  that 
for  circularly  symmetric  beam  distributions  the  radial  and  azimuthal  current  densities 
go  to  zero  as  r  — >  0,  the  Bessel  function  expansions  can  be  expressed  as 

p=jrpm{z,t)J0(i^\  (5.1a) 

m=  1  '  ' 

oo  /  .  \ 

(5.1b) 

m=  1  '  ' 

oo 

Jo  =  22 

m= 1 

and 

Jz  =  f2JZm(z,t)J0(J-^V  (5. Id) 

m=  1  '  ' 

where  Ji(x)  is  the  Ith  order  Bessel  functions  of  the  first  kind  and  jim  is  the  mth 
zero  of  Ji(x).  We  note  that  while  three  of  the  above  source  functions  utilize  jim  in 
connection  with  Ji(x)  as  Fourier- Bessel  series  expansions,  the  radial  current  density, 
Jr,  does  not.  In  particular,  Jr  utilizes  jom  zeroes  in  combination  with  J\  (x),  which  is 
a  special  type  of  expansion  known  as  a  Dini  series  [51]. 

The  coefficients  in  each  expansion,  which  are  functions  of  z  and  t,  can  be  im¬ 
mediately  found  using  orthogonality  relations  of  the  Bessel  functions  and  are  given 


Jimr 


(5.1c) 
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by 

Pm  =  2  T2(  ■ — 7  /  drrp(r,z,t)J0  (5.2a) 

a2Jf  (Jo m)  Jo  V  a  J 

J rm  =  ~TWr '■ — \  [  drrJr(r,z,t)J i  (5.2b) 

a?J(  (jom)  Jo  \  «  J 

Jem  =  2  J2(  ■ — 7  [  drrJe{r,z,t)Ji  ( )  ,  (5.2c) 

a2  Jo  (Jim)  Jo  \  a  J 

and 

Jzm  =  2  j2(  ■ — 7  [  drrJz(r,z,t)J0  (5. 2d) 

o2J\  \Jom)  Jo  \  «  / 

Since  Jgm  is  independent  of  6  coordinates,  the  continuity  equation  can  be  represented 
in  terms  of  pm,  Jrm,  and  Jzm,  i.e 


9 Pm 

~df 


jom  j  !  9 J zn 

'Jrm  i  7: 

a  az 


=  0. 


(5.3) 


In  addition,  Maxwell’s  equations  leads  to  inhomogeneous  wave  equations  for  the 
electromagnetic  holds  with  charge  and  current  densities  acting  as  drive  terms,  i.e., 

1  <92E  1  dJ  .  A  , 

V45  -  f  aF  =  7aVp  +  (5'4a) 

and 

1  <9^B 

v2b  -  f  aF  =  -»>v  x  J-  (54b) 

From  Eqn.  (5.4),  one  can  immediately  find  that  under  cylindrically  symmetric  con¬ 
ditions,  Jg  will  generate  the  fields  Eg ,  Br,  and  Bz,  while  p,  Jr,  and  Jz  generate  the 
other  holds  Er,  Ez,  and  Bg. 


5.1.2  TE  Mode  Space-Charge  Fields 

From  Eqn.  (5.4),  the  electric  and  magnetic  helds  generated  by  Jg  can  be  expressed 
in  terms  of  Bessel  functions  with  unknown  coefficients.  These  coefficients,  which 
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depend  on  z  and  t,  can  be  found  by  solving  wave  equations  with  Jg  as  a  source.  The 
expansions  of  TE-mode  fields  are  given  by 


Ee  =  Y,E0m{z,t)J1{Jlmr 

Br  =  J2Br™(z,t)J:  (jlmr 


m—  1 


and 


^  ^  t)  Jq 


m=  1 


Ji  mr 


(5.5a) 

(5.5b) 

(5.5c) 


At  r  =  a,  we  note  that  Eg  and  Br  equal  to  zero.  Eqns.  (5.5a)  and  (5.5b)  imply  that 
the  series  expansions  of  Eg  and  Br  automatically  satisfy  the  conducting  boundary 
conditions,  namely,  E\\ |  f  =  0  and  -Bj_|surface  =  0.  In  addition,  Eqn.  (5.5c)  gives 
another  boundary  condition,  dBn/drl  ,  =0. 

The  coefficients  in  the  held  expansions  are  derived  directly  from  Eqn.  (5.4)  and 
are  given  by, 


and 
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(5.6a) 

(5.6b) 

(5.6c) 


One  can  fold  the  solutions  of  the  Eqn.  (5.6)  using  time-dependent  Green’s  function 
method.  This  method  was  already  discussed  in  detail  for  longitudinal  currents  in  a 
conducting  pipe  in  Chapter  3.  In  particular,  the  solutions  to  the  held  coefficients  are 
given  by, 


Egm 


dtf  f  dz'/io  —Glm  (z,t;z',t') 


(5.7a) 
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Brm=  I  dt'  j  dz'no  dJem^  1  -  -  Glm  (z,  t]  z'l  t')  (5.7b) 

and 

Bzm  —  f  dt'  f  dz'no— j0m(z',t')Gim(z,t-z,,t>)  (5.7c) 

J  —  OO  J  ^ 

where 

Gim  (z,  t]  z',  f)  =  C-J0  ©  (A2)  (5.8a) 

is  a  time-dependent  Green’s  function  that  satisfies  the  differential  equation 

d2Gim  1  d2Gim  f  jlm\2  ^  x  (  !\XU  4->\  tr.  ok\ 

Glm  =  -S(z-z)S(t-t),  (5.8b) 

and  A  —  \j ' c1  (t  —  t')2  —  (z  —  z')2. 

Therefore,  for  a  given  Jg(r,z,t),  Jgm(z,t)  can  be  found  by  Eqn.  (5.2c),  and  then 
coefficients  Egm,  Brm,  and  Bzm  can  be  computed,  finally  the  TE  mode  electromagnetic 
fields  can  be  obtained. 


5.1.3  TM  Mode  Space-Charge  Fields 

We  now  turn  our  attention  to  the  TM  mode  fields,  Er,  Ez,  and  Bg.  In  a  similar  way, 
we  can  make  series  expansions  for  Er,  Ez,  and  Bg ,  from  Eqn.  (5.1)  and  (5.4).  i.e., 

OO 

Er  =  ^  ^  Errn(z,  £)  J\ 

m=l 

Ez  =  jr  Ezm(z,t)J0  ,  (5.9b) 

m= 1  '  ' 

and 

OO 

Bg  —  N  ^  Bgm(Zi  t)J\ 

m  1 

Again  note  that  the  held  expansion  for  Ez  correctly  satisfies  the  conductor  pipe 
boundary  condition  at  r  =  a,  namely,  E|||  =  0.  These  expansions  when  plugged 
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into  the  field  wave  equations,  i.e.,  Eqns.  (5.4a)  and  (5.4b),  yield  differential  equations 
for  the  field  coefficients  given  by: 


d2Erm  _  1  d2Erm 
d  z2  c 2  d  t2 


JO m  Pm 
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dt 
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8  Jr 


I1  O' 


dz 


(5.10a) 

(5.10b) 

(5.10c) 


In  a  similar  fashion,  we  can  construct  Er,  EZ1  and  Bg  using  a  second  time- 
dependent  Green’s  function.  In  particular,  the  coefficients  of  the  fields  are  given 
by, 

r1  r  (  l  n  n  i  \ 
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Erm 


7  j  /  7  /  I  Jo m  Pm  O'Jrm 

dt  I  dz  ( - /x0— — - 

a  6q  ov 


EZm 


Bf)rr>.  — 


dt'  /  dz 


,  /  1  dpm  .  dj, 

7a7F  +  m~ 


where 


dt' 
d  7 

7 1  /  /  7/1  JUm  T  .  uurrn 
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|  (*0m  (%)  t)  %  5  t  )  , 
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is  a  time-dependent  Green’s  function  that  satisfies  the  differential  equation 
d2Gnm  1  d2Gc\  x2 


Otn.  1  (7  Crom,  /  Jorr 

dz2  c 2  dt 2  V  a 


G0m  =  -5(z  -  z')  5  {t-  t') . 


(5.11a) 

(5.11b) 

(5.11c) 

(5.12a) 

(5.12b) 


Finally,  the  solutions  to  the  field  coefficients  in  Eqns.  (5.7)  and  (5.11)  along  with 
the  expansion  formulas  in  Eqns.  (5.2),  (5.5),  and  (5.9)  offer  a  complete  description 
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of  electromagnetic  space-charge  fields  for  cylindrically  symmetric  beams  in  a  circular 
conducting  pipe. 

We  should  note  that  for  the  expansions  given  in  Eqns.  (5.2),  (5.5),  and  (5.9),  there 
are  other  useful  formulas  connecting  the  field  coefficients  which  can  be  immediately 
found  from  Maxwell’s  equations.  In  particular,  one  can  easily  verify  the  following 
relations, 
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J  Om  J-y  _  T 
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1  dEzm 
c 2  dt 


(5.13f) 


One  can  immediately  see  that  although  all  of  the  field  coefficients  will  be  de¬ 
termined  from  integrations  of  Green’s  functions  over  z'  and  t'  for  given  charge  and 
current  densities  in  Eqns.  (5.7)  and  (5.11),  it  is  only  necessary  to  compute  these  inte¬ 
grations  for  a  few  of  the  field  coefficients.  Hence,  using  relations  in  Eqns.  (5.13),  it  is 
easy  to  determine  the  other  coefficients.  This  will  be  demonstrated  in  the  numerical 
example  in  the  next  Section. 
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5.2  Numerical  Implementation 


5.2.1  Description  of  the  System 

We  now  demonstrate  how  to  numerically  model  the  space-charge  fields  of  a  radially 
oscillating  bunched  beam  in  a  circular  pipe.  We  assume  that  the  beam  is  parabolic 
in  the  transverse  direction,  and  is  oscillating  in  the  radial  direction.  Specifically,  the 
beam  distribution  is  given  by 

p(r’  : 0 = t1  _  w>] 0  wt)  ■ r]  [e  {z + f )  ■ 0  (2  ~  f )  ]  ■  (5-14) 

where  Q  is  the  total  bunch  charge,  L  is  its  longitudinal  length.  The  beam  radius,  r&(t) 
is  a  constant,  Tq,  for  time  t  <  0,  and  then  is  sinusoidally  oscillating  with  amplitude 
Sr  and  angular  frequency  u  in  the  radial  direction  for  t  >  0,  i.e., 


r0,  for  t  <  0, 

r0  +  Sr  —  <5r  cos(cat),  for  t  >  0. 


(5.15) 


5.2.2  Derivations  of  the  Space-Charge  Fields 


For  given  charge  and  current  densities,  p(r,z,t )  and  J(r,z,t),  one  can  compute  the 
exact  space-charge  fields  Er,  Ez,  and  Be  using  the  Green’s  function  methods  discussed 
in  the  previous  section.  These  space-charge  fields  are  fully  self-consistent  and  include 
the  time-dependency  of  the  oscillations,  the  initial  transient  start-up  physics  at  time 
t  —  0,  and  the  effect  of  the  image  charges  and  image  currents  on  the  circular  pipe. 

By  specifying  the  charge  density  in  Eqn.  (5.14),  we  are  able  to  compute  the  charge 
density  coefficients,  pm,  in  Eqn.  (5.2a),  which  are  given  by 


rb/a,  drb/cdt 
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Figure  5.1:  Plots  of  and  ^  as  functions  of  time. 


5.2  Numerical  Implementation 


91 


Next,  we  calculated  the  coefficient  of  the  longitudinal  electric  held  Ezm  from 
Eqn.  (5.11b),  using  pm  in  Eqn.  (5.16).  Note  that  we  have  made  use  of  the  fact 
that  Jzm  =  0  since  no  longitudinal  currents  exist  when  implementing  Eqn.  (5.11b). 
From  Eqns.  (5.11b)  and  (5.12a),  Ezm  can  be  derived  as  following: 


Ezm  =  -jt  dt' J dz'  ^-dpm^yt')  Go™ CM; *',*') 

4 Q  [f  ( cdd_\  1  a2  r  2  / .  n\  _  / .  r&\ 

ne0a,L  J \  a  J  j%mJf(jom)  rf  [j0mrb/a  1  v  "  a  )  0  \J°m  a  J 

x  j„  (ifK)  e  (Ap  -  _)  e  (XI)  , 

(5.17) 


where  X±  =  \J c2  (t  —  t')2  —  (z  ±  L/2)2 .  Therefore,  the  longitudinal  electric  held  is 
given  by 

_  _  4 Q  x  ^  j'1  / cdE\  J0  (jomr / a)  a 2  f  2  /  .  ryx  _  /  .  7yv‘ 

7re0aL  ^  \  a  J  j$mJf  (j0m)  r\  [j0mrb/a  1  v0m  a  /  0  v0m  a  ) 

x  J0  f ifx+)  e  (Aq  -  j„  IA_)  e  (a2  )]  , 


(5.18) 


Since  remains  constant,  r0,  for  t  <  0,  we  can  manipulate  the  time  integration 
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in  the  following  manner: 


Ez  = 


4 Q  ST''  Jo  (jomr / o) 

^  J'L  Jl  (jom) 


r>0  cdt'  a 2  2 

-oo  a  rb  UomH/a 


I"  cdt' a2  2  /.  r0\  /.  r0\ 

/  9  /  Jl  I  JOm  I  Jo  l  JO m.  I 

/o  a  ro  Uomro/a  V  a/  V  a/J 

['  cdt'  a2  [  2  / 

—  / - 2  - - Jl  ( Jo 

Jo  a  ro  Uomro/a  V 


T  ( ■  rb\  j  (■  n\ 

J\  yjom.  q  J  Jq  yJOm  ^  J 

['  cdt'  a2  [  2  / .  r6\  /  .  r6\ 

/  2  /  "'ll  JOm  j  Jo  I  JOm  ) 

/o  o  rb  Uomn/a  V  aJ  \  aJ 


T  (  ■  r°\ 

Jo  yjom  a  J 


J„  J—K  e  (a2)  -  Jo  =a_  e  (a2  )  , 


(5.19) 

where  we  break  up  the  integration  into  the  two  terms  around  t  —  0,  and  add  and 
subtract  the  same  integration  term  with  the  radius  of  r0  instead  of  77, (t).  And  then 
we  join  the  second  and  fourth  terms  together,  and  do  the  same  for  the  first  and  third 
terms.  This  yields 


Er=- 


4 Q  ■sr-'  Jo  (jorn'C / a) 

7re°aX  m  ti0™) 

f*  cdt'  j  a2  2 
Jo  «  1  rl  [jornn/a 


T  ( ■  rb\  T  ( ■  rb\ 

J 1  yJOm  J  Jo  yJOm  ^  J 

a2  \  2  /.  r0\  /. 

2  ■  /  “1  l  JO in  )  Jo  l  JOn 

»o  LJom'ro/a  v  a/  V 

T  f  ■  r°\  T  (  ■  Ml 

Jl  I  Jo m  )  Jo  l  Jo m  ) 
i  \  a  J  \  a  J 


(5.20) 


['  cdt' a2  [  2  /.  r0\  /.  r0\l 

+  / - 9  - - 7  Jl  l  JOm —  ~  Jo  [JOm —  j 

J- 00  0  ro  LJom'ro/a  V  a/  V  a/J 

X  J„  (^A+)  0  (A2 )  -  Jo  (J==A_)  0  (A2,)  . 

This  method  is  useful  since  the  second  integration  term  is  exactly  equivalent  to  the 
electrostatic  part  of  Ez,  hence  there  needs  no  integration  up  to  a  time  t  =  —00.  The 
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field  will 


derivation  of  the  electrostatic 
term,  the  final  Ez  field  is  given  by 

4 Q  Jp  (jOmT  /  Q 


E. 
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for  z  >4: 


-*£“■(*+#) 

_30rn  ( 30  m  ( z  L\  ~|  ,  ,  r 

Jz+2i-e  2J  for  1 2 1  <  | 
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—  e~ 


22jnbz+L 


for  z  <  —  ^ 


(5.21) 


One  can  build  np  the  transverse  electromagnetic  field  components,  and  fig, 
using  Eqns.  (5.11a)  and  (5.11c)  with  pm  and  Jrm.  However,  calculations  of  the  time 
integrations  for  different  integrands  requires  large  computational  resources.  As  men¬ 
tioned  in  the  previous  section,  we  can  reduce  these  tasks  by  using  the  relations  from 
the  Maxwell’s  equations  (Eqn.  (5.13)).  For  given  charge  and  current  densities,  it  is 
easy  to  find  at  least  one  field  coefficient,  for  example,  Ezm,  in  this  numerical  example, 
as  well  as  those  of  charge  and  current  densities.  Then,  we  can  use  these  coefficients 
to  derive  other  components. 

Henceforth,  the  transverse  space-charge  electric  field  coefficients  Errn  and  Bgm  can 
be  found  from  Ezm  and  pm  using  Eqns.  (5.13c)  and  (5.13f),  respectively,  i.e., 


(5.22a) 
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a  [  t  |  f  dEzm 
Be  ~  2^  —  \ llJ zm  +  ^2  ~gf~ 


J\  (  jom 
\  a 


XA 


a  1  d  E~ 


j0m  c2  dt 


-/,  (.hJ-)  • 


(5.22b) 


Note  that  we  have  made  use  of  the  fact  that  Jzm  =  0  since  no  longitudinal  currents 
exist  when  implementing  Eqn.  (5.22b).  Using  the  relation, 


Jo  Of  aA  e  (\D 


3  T  ( jom  , 
YzJ°  [a 
d\±  d  r 

'  dz  d\±  [ 
_z  ±  L/2  'jo. 
A±  c 


Jo(^A±j0(A|) 


Ji  —  A±  ©  (A|)  -  2A±  Jo  3—\±  5  (A|) 


(5.23) 
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the  radial  electromagnetic  space-charge  held  is  derived  as 

Er 


8Q  J\  ( jomT / a) 

neoaL  0  m) 


4  Q 


a  i  (  .  rb 

rb  \_3omn/a  V  a 

cdt'  Ji  (j0mr/a) 


do  (  jom. 

a 


9 1 2  +  I 


eu-l 


E 


neoaL  a  jomdf(j0m) 


x 


2  I  '  /  ^1  ( Join 

rb  [jornn/a  V  a 


rb 


do  jo 


rb 


r  o 


x 


z  +  L/2 

z-L/2 

A^ 


di 


a 

Jo_ 

a 


Ji  (  4^A_  )  0  (A2 ) 
2 


2  I  '  /  ~d\  (  join 

Tq  l3omro/a  X  a 

2  a 


30m 
2  a 
JO  m. 


do  (jo 


r  o 


4Q  </] i  (join'C / a)  a2 


7re°aL  ^  JomJi(iom)  rjj  Liomro/a 


A+  )  0  (A2+)  -  ^A+ J0  A+ j  5  (A+) 

A_  J0  (^A_^)  5  (A?.) 

T  f  ■  ro\  j  ( ■  r0 

l  JO m  )  ^0  l  JO m 

\  a  /  V  a 


X  < 


_  JQm 


(  7A-  —  \  2Qrn(  ,_L_\ 

KZ+2  )  -|_  e  a  0  2  J 


•70m  ^  ~  |  ^  \  j()m  (  ~ E  A 

— e  a  2; 


for  z  >  -| 
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(5.24) 


In  the  derivation  of  Bg,  there  needs  more  delicate  consideration  since  both  in¬ 
tegrands  and  the  integration  range  in  the  time  integration  of  Ezm  depend  on  the 
differential  variable  t.  Let  the  integrands  be  /  (t,  t'),  then  the  time  integration  can  be 
simply  derived  as 

tt  f  dt'f  (t,  t')  —  f  (t,t)  +  f  dt'zrrf  (t,  t') .  (5.25) 

«/  — OO  */  — OO 

In  Eqn.  (5.25),  f(t,  t )  is  always  zero  since  A^  =  c2  (t  —  t')2  —  (z  ±  L/2 )2  are  less  than 
or  equal  to  zero  at  t'  =  t,  i.e.,  0  (A±)  do  not  contribute  to  the  calculation. 
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And  the  time  derivatives  of  the  Bessel  functions  in  the  integrands  are  similar  to 
Eqn.  (5.23)  as 

|Jo(^A±)e(Ai) 

4^Kta±)9®]  <5  26) 

=  _e^_T)  ^A±^J  0  (A|)  _  2A±J0  ^-X±^S  (A|)  . 

Therefore,  using  Eqns.  (5.9b),  (5.21),  and  (5.26) 

B  _  40  j'f  cdt'  Ji  ( jomr/a ) 

x(4Mv^-)-J .(^)1 

^6  |jom'r&/a  V  a/  V  a/J 

+  AT  (jo m~)  ~  Jo  (jom-)]  }  (5.27) 

Tq  [ Jom'r o/ a  V  a/  V  a/J  J  v  ' 

x  h  (if  A+)  8  M)  -  ^  (^)  *  M)] 

-  ^  [J‘  ( vA-) 8  h-)  -  (^-)  *  (A-)]  }■ 


5.2.3  Numerical  Requirements 

As  demonstrated  in  the  previous  chapter,  the  numerical  requirements  to  analyze  the 
electromagnetic  fields  are  studied.  The  number  of  radial  eigenmodes  needed  for  good 
convergence  scales  inversely  proportional  to  the  radial  size  of  the  perturbation  which 
is  being  analyzed.  Figure  5.2  shows  shows  normalized  Er  vs.  r/a  at  time  t  =  0.25T, 
where  T(=  27t/u ;)  is  the  characteristic  beam  oscillation  period,  i.e. ,  the  beam  edge  is 
at  r  =  ro  +  Sr  =  0.5a  for  ro  =  0.25a  and  Sr  =  0.25 a.  In  Figure  5.2,  the  normalized 
radial  electric  held  is  plotted  for  rn  —  5  radial  modes  (red:solid),  m  =  10  modes 
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Figure  5.2:  Plots  of  the  normalized  radial  electric  field  at  z  =  0 
and  t  =  0.25 T  are  plotted  using  5  radial  modes  (red:solid),  10 
modes  (greemdashed),  50  modes  (blue:dotted),  and  100  modes 
(pink: dashed- dotted)  (rmn  is  the  radial  eigenmode  number). 
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(green: dashed),  m  —  50  modes  (blue:dotted),  and  m  =  100  modes  (pink:dashed- 
dotted). 
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Figure  5.3:  Plots  of  the  percentage  error  of  the  normalized  radial 
electric  field  at  z  =  0  and  time  t  —  0. 


Figure  5.3  shows  the  percentage  difference  of  the  normalized  radial  electric  held 
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for  different  eigenmodes.  The  differences  are  calculated  from  m  =  2000  modes.  We 
note  that  as  we  increase  the  eigenmode  number  Er  is  rapidly  converged,  such  that 
the  difference  between  m  =  100  and  m  =  500  modes  is  already  negligible.  It  is 
obvious  that  most  discrepancies  are  occurred  at  r  =  0.0,  r  =  0.25a  =  r o  and  r  = 
0.5a  =  ry,  since  the  eigenmode  summations  are  conducted  around  these  radii  as  seen 
in  Eqn.  (5.24).  In  order  to  achieve  higher  accuracy,  we  use  m  =  1000  modes  for  the 
following  numerical  examples. 

In  Eqns.  (5.21),  (5.24),  and  (5.27),  numerical  time  integrations  needs  restrictions 
for  integration  time  step  size,  as  previously  discussed  in  Chapter  3.  The  time  step 
size  must  be  small  enough  to  resolve  the  radial  oscillation  of  the  beam,  which  is  rep¬ 
resented  by  the  Green’s  function  found  in  the  Eqns.  (5.21),  (5.24),  and  (5.27).  In 
Chapter  3  for  computing  the  electromagnetic  space-charge  fields  due  to  the  longitudi¬ 
nal  motion  of  the  beam,  we  found  that  the  time  step  size  is  restricted  by  the  relation 
At1  «  O.Olry/c.  In  order  to  achieve  highly  accurate  fields,  the  step  size  was  10”6  in 
normalized  units. 

When  we  compare  Eqn.  (3.38)  and  (5.24),  for  example,  the  temporal  part  in  the 
time  integration  is  very  similar,  although  the  beam  radius  in  the  present  work  also 
depends  on  the  integration  variable  t' .  Hence,  one  can  expect  that  the  required 
time  step  size  has  to  be  in  the  same  criteria.  However,  the  temporal  part  of  the  rf 
photoinjector  study  is  more  complicated,  i.e. ,  it  has  the  dependence  of  the  longitudinal 
beam  trajectory,  z"(t),  while  there  is  no  longitudinal  motion  in  this  model.  Therefore, 
the  requirements  of  the  time-step  size  for  the  integration  is  less  restricted  than  those 
in  the  previous  work.  Moreover,  for  the  parameters  in  this  example,  the  beam  radius 
is  varying  from  0.25a  to  0.75a,  so  that  a  time  step  of  At'  =  0.0001 a/c  is  enough  to 
get  excellent  numerical  convergence. 
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5.2.4  Space-Charge  Fields  with  Transverse  Current 

As  an  example  of  the  transverse  current  model,  we  show  the  plots  of  the  normalized 
radial  electric  held  vs.  r  at  the  middle  of  the  bunch,  z  —  0  and  the  normalized 
longitudinal  electric  held  vs.  z/a  at  the  axis,  r  =  0  and  vs.  r/a  at  z  =  0.001  for 
different  times,  which  are  related  to  the  characteristic  beam  oscillation  period,  T(= 
2tt/u).  We  choose  times,  t  —  0,  t  —  0.25 T,  t  =  0.5 T,  t  =  0.75T,  and  t  —  T,  in  which 
the  beam  is  oscillating  for  one  full  period.  The  plot  of  normalized  azimuthal  magnetic 
held  is  excluded,  since  it  is  zero  at  the  bunch  mid-plane,  z  =  0,  as  well  as,  at  the  axis, 
r  =  0,  clue  to  the  symmetry.  But  it  is  non-zero  at  different  r  and  z  values. 

In  order  to  hnd  the  electromagnetic  effects,  the  electrostatic  space-charge  holds 
are  also  calculated  and  are  plotted  for  the  comparison.  The  electrostatic  helds  are 
also  found  using  the  electrostatic  Green’s  function  [34],  and  the  detailed  derivations 
are  shown  in  the  previous  section.  The  beam  oscillation  starts  with  the  initial  beam 
radius,  =  0.25a,  the  radial  variation,  Sr  =  0.25 a,  and  the  oscillating  frequency  is 
u  =  joic/a,  where  a  is  the  pipe  radius  of  0.0908  m.  All  held  calculations  are  fulfilled 
with  the  numerical  requirements,  i.e.,  m  =  1000  modes  and  At'  =  O.OOOla/c. 

Figure  5.4  shows  plots  of  normalized  Er  vs.  r  for  different  times  during  one  char¬ 
acteristic  beam  oscillation  period.  The  red-solid  curve  is  the  electrostatic  held,  while 
the  blue-dashed  curve  is  the  electromagnetic  held.  Figure  5.4  has  plots  of  the  nor¬ 
malized  radial  electric  held  at  z  —  0.0  for  times  t  =  0  (Figure  5.4(a)),  t  =  0.25T 
(Figure  5.4(b)),  t  =  0.5 T  (Figure  5.4(c)),  t  =  0.75T  (Figure  5.4(d)),  and  t  =  T 
(Figure  5.4(e)).  Since  the  beam  is  not  oscillating  for  t  <  0,  the  two  electric  helds  at 
t  —  0  are  exactly  same,  i.e.,  fully  electrostatic,  as  shown  in  Figure  5.4(a).  At  time 
t  =  0.5T,  the  beam  has  the  maximum  radius,  ry  =  r'o  +  2 Sr  =  0.75a.  We  see  that  the 
solution  of  the  electrostatic  held  is  underestimated  by  about  10  %  compared  to  the 
electromagnetic  held  in  Figure  5.4(c).  For  t  —  T,  however,  one  can  see  that  the  elec- 
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(a)  t  =  0.0 


(d)  t  =  0.75T  (e)  t  =  T 


Figure  5.4:  Plots  of  the  normalized  radial  electric  field  at  z  =  0  and 
at  time  (a)  t  —  0,  (b)  t  —  0.25 T,  (c)  t  =  0.5 T,  (d)  t  =  0.75T,  and  (e) 
t  =  T  computed  electrostatically  (red:solid)  and  electromagnet ically 
(blue:  dashed). 
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trostatic  field  is  overestimated  compared  to  the  electromagnetic  field.  Since  the  beam 
distribution  is  parabolic  in  the  transverse  direction,  the  peak  is  located  at  slightly  less 
than  the  beam  edge.  When  we  compare  the  two  held  solutions  for  times  t  =  0.25T 
and  t  =  0.75T,  where  the  beam  radii  for  both  cases  are  rj,  =  tq  +  Sr  =  0.5a,  the 
electrostatic  held  is  underestimated  or  overestimated  around  the  peak.  However,  the 
peak  electromagnetic  held  for  t  =  0.25 T  is  less  than  the  peak  of  the  electrostatic  held, 
while  the  peak  electromagnetic  held  for  t  =  0.75T  is  larger  than  the  peak  electro¬ 
static  held.  This  is  because  the  beam  is  radially  oscillating  in  the  outward  direction 
for  t  =  0.25T,  while  it  is  moving  inward  for  t  =  0.75 T.  For  t  —  T,  one  can  see  that 
the  electrostatic  held  is  overestimated  compared  to  the  electromagnetic  held. 

Figure  5.5  shows  plots  of  a  normalized  Ez  vs.  z  at  r  —  0  (Figure  5.5(a))  and 
r  =  0.25 a  (Figure  5.5(b))  for  the  times  t  —  0,  t  —  0.25 T,  t  =  0.5 T,  t  =  0.75T, 
and  t  =  T.  From  Figure  5.5,  it  is  immediately  apparent  that  Ez  is  maximum  at 
t  —  0  and  t  =  T.  This  is  easily  explained  by  the  fact  that  as  the  bunch  radius 
increases,  the  beam  charge  density  at  r  —  0  decreases  leading  to  a  smaller  Ez.  We 
find  that  while  the  difference  in  Ez  as  found  by  electrostatic  and  electromagnetic 
methods  at  r  —  0  and  at  times  t  —  0  and  t  =  T  for  these  parameters  is  negligible, 
the  difference  in  Ez  is  approximately  20%  at  times  t  =  0.25T  and  t  =  0.75T  near  the 
longitudinal  beam  edges  and  increases  to  40%  at  time  t  =  0.5 T.  When  we  compare 
Figures  5.5(a)  and  5.5(b)  near  the  longitudinal  beam  edge,  the  electrostatic  solutions 
are  underestimated  on  the  beam  center,  r  —  0,  while  those  are  overestimated  as  r 
is  increased.  Therefore,  the  effects  of  the  longitudinal  electromagnetic  space-charge 
holds  are  important  near  the  beam  edge  for  the  radially  oscillating  beam  as  described 
in  this  example. 


5.2  Numerical  Implementation 


103 


z/a 

(b)  r  =  0.25 a 


Figure  5.5:  Plots  of  the  normalized  longitudinal  electric  field  at 
(a)  r  =  0  and  (b)  r  =  0.25 a  for  times  t  —  0,  t  —  0.25 T,  t  =  0.5 T, 
t  =  0.75 T,  and  t  =  T  computed  electrostatically  (red:solid)  and 
electromagnetically  (blue:dashed) . 
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Chapter  6 

Conclusions  and  Fixture  Plans 

In  this  dissertation,  we  have  presented  electromagnetic  space-charge  models  of  elec¬ 
tron  beams  in  the  presence  of  conducting  boundaries  in  accelerators.  We  utilized 
Green’s  function  methods  to  derive  analytical  solutions  of  the  electromagnetic  space- 
charge  potentials  and  fields  in  the  pipe  structures.  We  have  developed  a  novel  com¬ 
putational  code,  called  IRPSS,  to  compute  these  fields  numerically.  Specifically,  we 
calculated  the  space-charge  fields  in  the  RF  photoinjector,  and  then  simulated  the 
beam  dynamics  near  the  cathode.  The  numerical  results  were  compared  with  those  of 
an  electrostatic  code.  We  found  significant  differences  between  electromagnetic  and 
electrostatic  space-charge  fields.  In  addition,  we  also  extended  our  Green’s  function 
techniques  by  including  transverse  currents.  We  applied  this  method  to  investigate 
the  electromagnetic  space-charge  fields  for  a  radially  “breathing  beam”  oscillation. 
At  the  end  of  this  chapter,  we  will  provide  future  plans  to  enhance  the  IRPSS  code 
for  modeling  the  rf  photoinjector  and  high-power  microwave  sources. 
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6.1  Space-Charge  Modeling  of  the  RF  Photoinjec¬ 
tor 

In  Chapter  3,  we  have  presented  analytical  solutions  of  the  electromagnetic  fields 
inside  of  a  pipe  structure  with  the  flat  cathode  at  z  =  0  using  a  compact  Green’s 
function  method.  These  solutions  include  effects  of  image  charges  and  currents  due 
to  the  cathode  boundary  and  the  necessary  causality  condition  due  to  the  finite  speed 
of  light.  The  accuracy  of  this  method  is  verified  using  a  benchmark  study. 

As  an  example  to  model  the  rf  photoinjector,  we  showed  the  exact  solutions  of  the 
electromagnetic  space-charge  fields  within  the  cylindrical  pipe.  The  code,  IRPSS,  has 
an  ability  to  calculate  the  space-charge  holds  in  the  rf  photoinjector.  We  found  the 
necessary  numerical  requirements  for  eigenmode  summations  and  time  integrations 
to  achieve  high  accuracies,  i.e.,  <  1  %  errors.  To  model  a  finite-size  bunched  beam, 
a  multi-slice  method  was  introduced. 

In  order  to  resolve  the  electromagnetic  fields,  the  distribution  function  of  the  zero¬ 
thickness  beam  was  expanded  in  terms  of  transverse  eigenfunctions.  We  estimated 
the  required  eigenmode  number  for  a  given  charge  distribution.  The  minimum  eigen¬ 
mode  number  depends  on  the  radial  beam  size  of  the  system.  For  example,  at  least 
2000  eigenmodes  are  needed  to  get  accurate  results  for  the  BNL  2.856  GHz  rf  pho¬ 
toinjector.  The  time  integrations  are  conducted  with  the  time  evolution  of  the  beam. 
The  numerical  simulation  showed  that  the  time  integration  step,  cAt'/a,  must  be  also 
much  less  than  a/c. 

Since  the  real  beam  has  a  finite  longitudinal  bunch  length,  the  multi-sliced  bunch 
model  is  introduced.  Each  slice  with  zero-t hickness  has  an  equal  amount  of  charge, 
but  total  charge  is  conserved.  We  found  the  required  number  of  slices  to  express  the 
finite  size  beam.  For  a  9  ps  bunch  length  withing  the  AWA  gun,  we  modeled  the 
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longitudinally  uniform  beam  with  the  multi-sliced  beam,  using  41  slices. 

6.2  Application  to  the  RF  Photoinjector 

The  computational  code,  IRPSS,  was  developed  to  solve  the  electromagnetic  space- 
charge  fields  and  to  track  the  electron  beam  in  the  rf  photoinjector.  First,  we  exam¬ 
ined  how  the  electromagnetic  space-charge  fields  differed  from  the  electrostatic  fields 
near  the  cathode. 

The  electrostatic  space-charge  fields  were  calculated  in  two  ways:  (1)  the  dis¬ 
cretized  bunch  method  with  multi-slices,  which  is  similar  to  the  electromagnetic  field 
calculations,  and  (2)  the  uniformly  finite  bunch  with  constant  Jz.  These  two  methods 
agreed  very  well,  and  verified  our  electrostatic  results.  The  comparisons  between  elec¬ 
tromagnetic  and  electrostatic  forces  showed  that  there  are  significant  discrepancies  in 
transverse  space-charge  force  near  the  cathode,  and  the  differences  were  decreased  as 
the  beam  was  moving  away  from  the  cathode. 

In  addition,  we  studied  the  effect  of  time-retardation  in  the  rf  photoinjector  during 
the  emission  of  the  electron  beam.  This  study  was  done  by  measuring  the  bearnloss 
both  numerically  with  electromagnetic  and  electrostatic  space-charge  field  solvers 
and  experimentally  in  the  ANL  AWA  1.3  GHz  rf  photoinjector.  Discrepancies  were 
found  between  two  critical  space-charge  fields,  especially,  for  beams  with  longer  bunch 
length. 

6.3  Transverse  Current  Modeling 

Chapter  5  extended  the  previous  Green’s  function  technique  including  the  effects  of 
the  transverse  beam  currents.  We  assumed  that  the  beam  is  cylindrically  symmetric 
and  the  geometry  of  the  conducting  cavity  system  is  also  cylindrically  symmetric. 
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With  the  help  of  these  symmetries,  we  expanded  the  beam  charge  and  current  densi¬ 
ties,  and  the  electromagnetic  space-charge  fields  in  terms  of  series  of  Bessel  functions. 
These  formed  a  set  of  differential  equations  from  Maxwell’s  equations,  and  we  found 
the  solutions  using  the  time-dependent  Green’s  function  methods.  The  Green’s  func¬ 
tions  automatically  satisfied  the  boundary  conditions. 

Using  this  formalism,  we  derived  the  analytical  solutions  of  the  electromagnetic 
space-charge  fields  in  the  presence  of  azimuthal  beam  currents  or  radial  and  longitu¬ 
dinal  beam  currents.  As  a  numerical  example,  we  computed  the  space-charge  fields 
for  a  radially  “breathing”  bunched  beam,  which  can  be  realizable  in  a  klystron  drift 
tube.  When  we  compared  the  electromagnetic  fields  with  those  of  the  electrostatic 
method,  the  differences  between  these  two  fields  were  non-negligible.  The  main  dif¬ 
ferences  of  the  longitudinal  space-charge  fields  occurred  near  the  beam  edge  up  to 
20%.  Moreover,  the  differences  in  the  transverse  space-charge  fields  were  varying  as 
the  beam  is  radially  moving  inward  or  outward.  These  numerical  comparisons  tells  us 
that  our  space-charge  formalism  may  play  an  important  role  in  modeling  space-charge 
dominated  beams  in  high-power  microwave  sources,  such  as  a  klystron. 

6.4  Future  Plans 

In  the  modeling  of  rf  photoinjectors,  we  simplified  the  geometry  of  the  system  only 
with  the  conducting  cavity  and  the  flat  cathode.  We  are  planning  to  extend  the 
IRPSS  code  by  including  the  effects  of  iris(es)  or  discontinuities  of  the  cavity.  These 
effects  can  be  applied  by  various  approaches,  however,  we  are  currently  focused  on  the 
method  developed  by  H.  Bethe.  This  method  includes  the  effects  of  the  hole  between 
two  consecutive  cavities  using  the  perturbative  approximations  of  field  expansions. 
To  the  lowest  order,  the  hole  effects  are  assumed  to  be  equivalent  to  oscillating  electric 
and  magnetic  multipole  moments.  The  superposition  of  the  field  expansion,  which 
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is  computed  without  iris(es),  and  the  multipole  moments  will  represent  the  total 
electromagnetic  space-charge  fields. 

The  space-charge  field  calculations  with  the  IRPSS  code  in  the  current  stage  are 
less  self-consistent  since  we  used  the  prescribed  beam  trajectory,  which  is  calculated 
in  the  presence  of  the  external  rf  field.  This  assumption  is  valid  when  we  examine  the 
beam  dynamics  near  the  cathode.  However,  the  beam  distribution  will  be  changed 
as  the  beam  passes  through  the  gun.  We  will  improve  the  code  to  self-consistently 
calculate  the  trajectories  due  to  both  the  external  fields  and  the  space-charge  fields. 

We  are  also  planning  to  study  how  the  beam  dynamics  can  be  affected  due  to  the 
electromagnetic  space-charge  fields  in  the  designs  of  magnetic  focusing  schemes  for 
photoinjector  guns,  such  as  emittance  compensation  methods.  This  study  will  inves¬ 
tigate  how  the  fields,  such  as  rf,  solenoid,  and  space-charge  fields,  are  coupled.  More 
improved  experimental  measurements  will  be  performed  to  see  the  time-retardation 
effects  in  the  rf  photoinjector. 

In  addition,  we  are  developing  space-charge  field  solvers,  which  can  include  arbi¬ 
trary  beam  currents,  for  example,  azimuthally  varying  currents.  Moreover,  we  are 
also  developing  the  analytical  techniques  for  including  more  complicated  conduct¬ 
ing  geometries,  such  as  coaxial  structures,  and  for  various  beam  charge  and  current 
distributions,  such  as  annular  beams. 
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Appendix  A 

Electromagnetic  Space-Charge 
Fields 


In  this  chapter,  we  will  derive  the  electromagnetic  space-charge  potentials  and  fields 
using  Green’s  function  methods  for  the  circularly  symmetric  pipe  with  the  flat  cath¬ 
ode.  We  assume  that  the  beam  has  only  longitudinal  currents.  Then,  the  charge  and 
current  densities  of  the  zero-thickness  beam  are  given  by 

p{r,t)  =  ^  (l  -O  ©  {rb  —  r)  S[z  —  z"(t)\ ,  (A.la) 

7T  T'b  \  rb  J 

and 

Jz  (r,t)  =  ^  (l  -  ^  0  {rh  -  r)  S[z-  z"(t )]  (L  ^  ,  (A. lb) 

7Trb  \  rb  J  at 

where  Q  is  the  bunch  charge,  is  the  beam  radius,  and  z"{t)  is  the  slice  location  at 
time  t.  We  normalize  the  variables  by  the  pipe  radius,  a,  and  the  speed  of  light,  c, 
such  as,  f  —  r/a,  z  —  z/a,  and  r  =  ct/a.  All  normalized  quantities  are  dimensionless. 
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A.l  Scalar  and  Vector  Potentials 


The  scalar  and  vector  potentials  can  be  obtained  by  the  Green’s  function  methods: 

0  (r ,t)  —  —  f  I  G<j)  (r,  t;  r',  t')  p  (V,  t')  d3r'dt',  (A. 2a) 

^0  J — oo  J 

and 

rt  r 

(A. 2b) 


Az  (r,  t)  =  n o  /  I  Ga  (r,  t;  r',  t')  Jz  (V,  t')  d3r'dt', 

where 

G<t>  =  %'52'l>mn{r±)ll>mn{r,±) 


Jo  ©  (Aq_)  -  Jo  ©  (A§+) 


ga  =  t:Y1  ^mn  (r±) 


A0- 


Jo  (  j  ©  (Aq_)  +  Jo  ©  (Ajj+) 


(A. 3a) 


ipmn  (r±)  = 


.imO 


1  Jm  {jmnX  /  J)  & 

|  Jm+1  (. j  ran)  | 


and 


/\2 


Aq±  =  c2  (t  —  t')  —  (z  ±  z') 


(A. 3b) 
(A. 3c) 

(A. 3d) 


With  the  simple  algebraic  calculation,  the  scalar  and  vector  potentials  are  given  by 


JJ  ^  4  cQ  ^  Jo  (jon  ra ) 

0(r,  z,  t)  = - 2  V  -2  I  T  7  ■  m 

rl  “  JonlJlUOn)r 


.jonTb/ O' 


J\  (jon~~)  ~  J°  (Ao 


n=0 

X  A_\  _  2  /  A+  i  ,  2 


'  T± 

a 


x  /  d,t' 

J  —oo 


Jo  (  JW—  J  ©(^_)  —  J0  (  J°n~  )  ©(^+) 


4Q  _ Jo  (jorj  ) 


-  Jl  (jorAb)  -  Jo  (jOrAb) 


^O'^O  ^  jon\JlUon)\2  [jonh 

Jo  (jon^-)  0(A0)  —  Jo  fjonA+)  ©(A+) 


(A. 4a) 


x  /  dr 

J  — OO 
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and 


4  (  ±\  ^cQdo  Jo  (jo nra) 

Az(r,z,t)=  __2 


*rb  ^  ion  I  J1  (ion)  | 2  Lion^ 


.  Tb 


J\  (ion —  )  —  Jo  (i On  — 


.  U 


x  /  df' 


A_ 


do  (  ion -  ©(A_)  +  Jq  I  jon -  ©(A.) 


Aj 


dz"  ( t ) 
dt 


J  t=t' 


^cQuo  Jo(jonf) 


AT 


^6°  ^ionlJl(ion)|2  Lion^ft 


Jl  (jonTb)  Jq  (jon.j’b) 


(A. 4b) 


x  /  dr' 


Jo  ( ion  A—  )  @(A2)  +  Jq  (ion.A+ )  ©(Ai) 


dz" (r) 
dr 


where  A^_  and  A^  are  given  by 


A2  =  c2  (t  —  t')2  —  [z  ±  z"(t'jf  ,  (A. 5a) 

and 

A ±  =  (t  —  r')2  —  [i  ±  A"(r')]2  .  (A. 5b) 


The  electromagnetic  helds  can  be  obtained  from  these  potentials. 


E  =  -V0  - 
B  =  V  x  A. 


dA 
dt  ’ 


(A.6) 
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Using  Eqns.  (A. 4)  and  (A.6),  it  is  easy  to  derive  the  transverse  electric  and  magnetic 
helds,  and  they  are  given  by 


Er(r. ,  z,  t)  •  ^ 

4  Q 


J\  (jonf) 


y  _ 

7T eor^a2  y  j0n\ Jl(ion)|2  IjOnU, 


Jl  (jonTb)  -  Jo  (jonh) 


(A.  7a) 


x  /  dr' 

J  —  oo 


Jo  (ionA-  )  0(A_)  —  Jq  (  JOnA+  )  0(A+) 
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Ft  A 

Bo(r,z,t)  = 


_  4cQ/i0  ^  -J\  (jo nr)  2  . .  ^  „ 

—  - ~ - - rrr - 7l2  - - “A  U0nr b)  —  J 0  (jOrA ft) 

JOn\MjOn)\2  [jOnH 
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(A. 7b) 


2  \1  r dz"(r) 


Using  the  relation,  c2  =  (/x0e0)  * , 


rw  JiUonf)  I”  2 

c  *r’*’  ’  ~  iceory  h  JMjon)]2  [jonh 


~J\  (jOn'U)  Jo  (jOrcTfe) 
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2x1  f  dz" (t) 


(A.8) 
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A. 3.1  Derivation  of  E~ 


The  longitudinal  electric  held  is  evaluated  from 


p  _  deft  d Az 

z~  dz  dt 


(A.9) 


And  we  will  derive  each  term  separately. 


T 

1  dz 


Since  the  time  integration  variable  t'  is  independent  of  z,  it  is  easy  to  take  a  partial 


derivative  of  eft  with  respect  to  z,  i.e., 
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Using  the  chain  rule  and  the  relation,  dJ0/dx  =  —J\(x), 
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For  this  calculation,  one  needs  to  be  careful  since  both  the  range  of  integration  and 


integrands  depend  on  the  variable  t. 

dAz  _  4cQ/i0  '  Jo  (jon~)  2 
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It  is  helpful  to  use  the  relation, 
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Using  the  chain  rule  and  the  property  of  the  Bessel  function, 
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Therefore, 
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A. 3. 2  Analytic  form  of  Ez 

Finally,  using  Eqns.  (A.  13)  and  (A.  19),  the  longitudinal  electric  field  is  given  by 
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Appendix  B 

Electrostatic  Space-Charge  Fields 


B.l  Electrostatic  Space-Charge  Fields  with  Delta- 
Function  Source 


Using  the  Green’s  function,  the  electrostatic  space-charge  potential  of  the  real  charge 
disk  in  the  beam  rest  frame  can  be  calculated  as  (See  pl43  in  Jackson): 

<j>(r,z)  =  —  [  dPr'G^  (r;  r')  p  (r') ,  (B.l) 
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(B.2a) 


(B.2b) 


and 

P  (r)  =  — \ 0  (rb  ~  r)  ( 1  -  r— )  J  (z  -  z0)  ■  (B.2c) 

nrb  V  rbJ 

As  seen  in  p  (r),  a  disk  beam  is  located  at  z  =  zo,  and  the  beam  motion  is  stationary. 
Since  the  image  disk  is  located  at  z  —  —  zq,  the  scalar  potential  is  given  by  the 
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superposition  of  potentials  of  real  and  image  disks: 
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Since  the  beam  is  at  rest  in  this  frame,  the  vector  potential  Az  is  zero. 

The  scalar  and  vector  potentials  in  the  laboratory  frame  can  be  calculated  with 
the  Lorentz  transformations,  as 

</>'  =  7  (0  +  c/3  •  A) , 


/ ^ 


(B.4a) 

A'  =7  [Az  +  ZjJ  ,  (B.4b) 

where  unprimed  potentials  are  in  the  rest  frame,  and  primed  potentials  are  in  the  lab 
frame.  From  these  potentials,  the  electrostatic  space-charge  fields  in  the  accelerating 
frame  can  be  derived,  i.e., 
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The  7  factor  in  front  of  z  is  for  the  Lorentz  transformation  of  ^-direction.  In  the 
B'e  calculation,  the  image  disk  is  moving  in  the  opposite  direction,  hence,  its  f3  has 
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the  opposite  sign.  We  assumed  that  the  beam  is  moving  only  in  the  ^-direction,  i.e., 
d/3/drj_  =  0.  Hence,  <9y/ dr  =  0,  but  dj/dz  ^  0.  The  Lorentz  transformation  of  the 
fields  (See  p558  in  Jackson)  are 
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and  Ez  is  not  changed  during  the  transformation: 
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B.2  Electrostatic  Fields  with  a  Constant  Jz 


In  this  section,  we  assume  that  the  bunch  has  a  finite  size,  i.e.,  it  is  not  broken  up 
into  a  number  of  slices,  and  the  current  density  is  constant  through  the  bunch.  This 
is  helpful  for  forming  an  analytical  electrostatic  result.  Then,  the  charge  density  is 
given  by 
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where  Q  is  the  total  bunch  charge,  r&  is  the  beam  radius,  and  z\  and  z 2  are  the 
locations  of  the  tail  and  head  of  the  bunch  at  a  specific  time,  respectively.  Using 
Eqns.  (B.l),  (B.2),  and  (B.9),  the  scalar  potential  is  computed  as 
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for  zi  =  0  and  Z2  =  Zb-  In  Eqn.  (B.10),  a  is  the  cavity  radius,  Erf(z)  is  the  error 
function,  and  Erfi  (2)  are  the  imaginary  error  function.  The  error  functions  are  defined 
by 

2  fz 

Erf(z)  =  —j=  /  dte~l  ,  (B.lla) 

v  n  Jo 

and 

Erfi(z)  =  —iErf(iz),  (B.llb) 


where  %  =  \f—l. 

In  order  to  include  the  image  charge  effect  due  to  the  cathode,  the  potential  with 
the  additional  bunch  with  the  opposite  charge  in  the  negative  z-direction  will  be 
added  such  that, 

0total  =  0real  T  0image-  (B.12) 


Since  the  longitudinal  electrostatic  space-charge  field  is  invariant  during  the  Lorentz 
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transformation,  El  =  EZ)Vest  =  —dcp/dz  in  the  laboratory  frame  and  is  given  by 
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nite  size  bunch  length 


In  this  section,  the  bunch  has  a  finite  length  in  the  longitudinal  direction  as 
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Then,  the  electrostatic  potential  is  given  by: 
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From  the  electrostatic  potential, 
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of  Korea,  12/2004  -  12/2006 


Professional  Experience 


•  Indiana  University  Cyclotron  Facility,  Research  Associate,  01/2005  -  Present 

•  Indiana  University,  Bloomington,  Associate  Instructor,  08/2007  -  12/2007 

•  Korea  Institute  of  Radiological  And  Medical  Sciences  (KIR AMS),  Research  As¬ 
sistant,  12/2002  -  12/2004 


USPAS  Participations 


Year  Courses  (Instructors) 

2001/Winter  Accelerator  Fundamentals  (undergraduate) 

(Michael  Syphers) 

2002/Winter  Accelerator  Physics 
(Gerald  Dugan) 

2003/Summer  Accelerator  Physics 
(SY.  Lee) 

2004/Winter  Classical  Mechanics  and  EM  for  Accelerators  and  Beams 
(Helmut  Wiedemann) 

2004/Summer  The  SNS  -  I,  Front  End  and  Linac 
(Thomas  Wangler  and  James  Billcn) 

2004/Summer  The  SNS  -  II,  Ring  and  Transport  Systems 
(Jie  Wei  and  Yannis  Papaphilippou) 

2005/Winter  Advanced  Beam  Dynamics 

(Bruce  Carlsten  and  Steve  Russell) 

2005/Winter  Computational  Methods  in  Beam  Dynamics 
(Robert  Ryne) 

2005/Summcr  Computational  Methods  in  Electromagnetism 
(Kwok  Ko  and  Zenghai  Li) 

2006/Winter  High  Current  Beam  Physics  in  Induction  Accelerators 
(George  Caporaso  and  Yu-Jiuan  Chen) 

2006/Winter  Fundamentals  of  Ion  Sources 

(Ian  Brown  and  Xavier  Godechot) 

2006/Summer  Beam  Physics  with  Intense  Space  Charge 
(John  Barnard  and  Steven  Lund) 

2007/Winter  Collective  Effects  and  Wakefields 
(Roger  Jones) 

2007/Winter  Modern  Computational  Accelerator  Physics 

(Panagiotis  Spentzouris  and  James  Amundson) 

2007/Summcr  Beam  Experiments  and  Measurements  at  the  NSCL 

(Marc  Doleans,  Guillaume  Machicoane,  Walter  Hartung,  and  Eduard  Pozdeyev) 
2007/Summer  Fundamentals  of  Low-Beta  Linear  Accelerators  with  Simulation  Lab 
(John  Staples  and  George  Gillespie) 

2008/Winter  Synchrotron  Radiation  and  Free  Electron  Lasers 
(Kwang-Je  Kim,  Yin-e  Sun,  and  Zhirong  Huang) 

2009/Winter  Special  Topics  in  Accelerator  Physics 
(Alex  Chao  and  Yunhai  Cai) 


Publications  and  Presentations 


•  C.S.  Park  and  M.  Hess  “The  Effects  of  Electromagnetic  Space-Charge  Fields 
in  RF  Photocathode  Guns,”  submitted  to  Nucl.  Instr.  Meth.  A. 

•  M.  Hess  and  C.S.  Park,  “Exact  Space-Charge  Field  Solutions  for  Beam  Cur¬ 
rents  in  a  Circular  Conductor  Pipe,”  submitted  to  Phys.  Rev.  ST  Accel. 
Beams. 

•  C.S.  Park  and  M.  Hess,  “Time-Retardation  Effect  Causing  Beamloss  in  the 
RF  Photoinjector,”  in  PAC’09  proceedings. 

•  M.  Hess  and  C.S.  Park,  “A  Dispersion  Free  Three-Dimensional  Space-Charge 
Modeling  Method,”  in  PAC’09  proceedings. 

•  C.S.  Park  and  M.  Hess,  “Studies  of  Electromagnetic  Space-charge  Fields  in 
RF  Photocathode  Guns,”  in  EPAC’08  proceedings. 

•  S.  Y.  Lee,  J.  Kolski,  Z.  Liu,  X.  Pang,  C.S.  Park,  W.  Tam,  and  F.  Wang, 
“Low  energy  electron  storage  ring  with  tunable  compaction  factor,”  Rev.  Sci. 
Instrum.  78,  075107  (2007) 

•  M.  Hess  and  C.S.  Park,  “A  Multislice  Approach  for  Electromagnetic  Green’s 
Function  Based  Beam  Simulations,”  in  PAC’07  proceedings. 

•  C.S.  Park  and  M.  Hess,  “Computational  Requirements  for  Green’s  Function 
Based  Photocathode  Source  Simulations,”  in  PAC’07  proceedings. 

•  M.  Hess,  C.S.  Park  and  D.  Bolton,  “A  Green’s  Function  Based  Space-Charge 
Field  Solver  for  Electron  Source  Simulations,”  Phys.  Rev.  ST  Accel.  Beams 
10,  054201  (2007) 

•  M.  Hess  and  C.S.  Park,  L.  Zhu,  and  D.  Bolton,  “Electromagnetic  Green’s 
Function  Based  Simulations  of  Photocathode  Sources,”  in  LINAC’06  proceed¬ 
ings. 

•  M.  Hess  and  C.S.  Park,  “Advanced  Electromagnetic  Analysis  of  Electron 
Source  Geometries,”  in  PAC’05  proceedings. 

•  J.-S.  Chai,  Y.-S.  Kim,  T.K.  Yang,  D.H.  An,  I.S.  Jung,  H.S.  Chang,  G.B.  Kim, 
J.S.  Kang,  C.S.  Park,  M.G.  Hur,  B.H.  Hong,  S.S.  Hong,  S.W.  Kim,  M.Y.  Lee, 
W.T.  Hwang,  J.H.  Kim,  H.H.  Shim,  W.Y.  So,  M.  Sekiguchi,  “Development  and 
Status  Report  of  Medical  Cyclotron  at  KIRAMS,”  in  Cyclotrons’04  proceedings. 


•  C.S.  Park,  J.S.  Kang,  D.H.  An,  J.S.  Chai,  M.  Sekiguchi,  “Beam  Dynam¬ 
ics  Study  of  High  Intensity  Compact  Cyclotron  for  Neutron  Generator  at  KI- 
RAMS,1'  in  Cyclotrons’04  proceedings. 


